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Abstract 

The  purpose  of  this  research  was  to  develop  a  method  for  sensor  placement  on  a 
Border  Patrol  interdiction  network.  Specifically,  this  thesis  sought  to  develop  a  proof  of 
concept  model  using  Microsoft  Excel,  with  some  add-on  capabilities,  to  optimize  the 
probability  of  detecting  intruders  who  have  already  breached  the  border  through  the 
placement  of  electronic  sensors  on  a  network.  A  model  was  developed  which  maximizes 
the  probability  of  detecting  intruders  by  optimizing  the  build-up  of  a  distributed  sensor 
network  subject  to  a  budgetary  constraint.  Several  different  optimization  algorithms  were 
developed  for  use  with  the  model.  All  were  tested  and  their  results  were  analyzed 
revealing  two  very  promising  sensor  placement  methods  for  optimizing  sensor  coverage 
on  a  network. 

Due  to  its  ease  of  use  and  ability  to  run  in  Microsoft  Excel,  it  is  believed  that  the 
model  developed  in  this  research  can  also  be  used  in  a  number  of  military  applications 
where  border  security  is  necessary. 
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OPTIMIZING  DISTRIBUTED  SENSOR  PLACEMENT  FOR  BORDER  PATROL 
INTERDICTION  USING  MICROSOFT  EXCEL 

I.  Introduction 

Background 

The  United  States  of  America  is  not  only  the  world’s  sole  true  superpower,  but  also  its 

most  hospitable  host.  Founded  by  immigrants  for  immigrants,  the  U.S.  takes  in  more  of  the 

world’s  poor  and  downtrodden  than  any  other  nation.  In  fact,  as  of  2006  “the  United  States 

accepted]  more  legal  immigrants  as  permanent  residents  than  the  rest  of  the  world  combined” 

(9).  In  America,  a  poor  immigrant  can  become  a  CEO  of  a  major  corporation  or  the  governor  of 

the  most  populous  state  in  the  nation.  A  study  from  Duke  University  found  that 

25  percent  of  technology  and  engineering  companies  started  from  1995  to  2005  had  at 
least  one  senior  executive  -  a  founder,  chief  executive,  president  or  chief  technology 
officer  -  born  outside  the  United  States.  (10) 

While  legal  immigration  has  always  been  a  great  boon  to  the  United  States,  recently  there 
has  been  an  alarming  increase  in  illegal  immigration.  The  government  of  the  United  States  has 
been  unable  or  unwilling  to  stop  the  flow  of  smugglers  and  illegal  aliens  across  its  borders.  Over 
the  years,  millions  of  people  have  entered  the  country  illegally;  mostly  across  the  southern 
border.  In  fact,  as  of  2004,  there  were  an  estimated  12  to  20  million  illegal  aliens  in  the  United 
States.  (11) 

In  the  past,  the  internal  debate  in  the  United  States  for  and  against  illegal  immigration  has 
centered  mostly  on  economics.  Businesses  have  enjoyed  the  cheap  labor  provided  by  illegal 
immigrants  while  workers  and  union  groups  have  decried  the  mass  hiring  of  illegal  aliens 
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(especially  in  sectors  such  as  agriculture  and  construction)  as  a  cheaper  alternative  to  American 
citizens  and  legal  immigrants. 

However,  with  the  rise  of  the  radical  Islamic  movement  in  the  1980s  and  90s,  and  the 

subsequent  terrorist  attacks  of  September  11,  2001,  overlooking  illegal  border  crossings  is  no 

longer  an  option.  Nowadays,  it  is  not  just  drug  smugglers  and  poor  people  looking  for  work  that 

are  of  concern.  There  is  now  a  very  real  threat  of  terrorists  using  our  porous  borders  to  infiltrate 

and  attack  our  country.  The  nation  must  protect  its  borders  in  order  to  protect  its  citizens  from 

the  very  real  threats  which  face  them.  For  this  reason,  the  U.S.  government  must  take  a  careful 

look  at  all  people  that  are  entering  this  country  and  do  its  utmost  to  prevent  unauthorized  entry 

into  the  United  States.  However,  despite  the  ever-present  threat,  it  is  estimated  that  an  average  of 

half  a  million  people  cross  the  U.S. -Mexico  border  illegally  each  year  and  the  U.S  Border  Patrol 

catches  only  1  in  4  border  crossers.  (5;  12:14) 

Recently,  in  an  effort  to,  among  other  things,  improve  the  management  of  its  borders,  the 

U.S.  government  has  undertaken  a  large  reorganization  effort.  On  March  1,  2003,  as  part  of  a 

realignment  effort  after  September  11,  2001, 

the  Department  of  Homeland  Security  was  established.  It  was  and  is  the  largest 
reorganization  of  our  Federal  Government  in  over  50  years.  As  part  of  the  Department  of 
Homeland  Security  reorganization,  U.  S.  Customs  and  Border  Protection  -  CBP  -  was 
created  by  unifying  all  frontline  personnel  and  functions  with  law  enforcement 
responsibilities  at  our  nation’s  borders,  that  is,  at  all  300  plus  ports  of  entry  of  the  United 
States  -  land,  sea  and  air  -  and  the  areas  in  between  the  official  ports  of  entry.  (13:2) 

After  September  11,  2001,  the  new  priority  mission  of  the  Border  Patrol  became  to 

“prevent  terrorist  and  terrorist  weapons  from  entering  the  United  States.”  Of  course  the 

traditional  mission  of  the  Border  Patrol,  to  prevent  “illegal  aliens,  smugglers,  narcotics,  and  other 

contraband,  from  entering  the  United  States,”  remains,  and  in  fact  coincides  with  this  new 

priority.  (13:2) 


2 


Unfortunately,  creating  a  new  bureaucracy,  with  central  leadership,  while  a  commendable 

idea,  does  little  to  actually  improve  the  desperate  situation  at  our  borders.  For  example,  in 

January  of  2005,  while  pursuing  three  suspected  drug  smugglers  in  SUVs, 

three  Hudspeth  County  [Texas]  deputies  and  at  least  two  Texas  Department  of  Public 
Safety  troopers  squared  off  against  at  least  10  heavily  armed  men  from  the  Mexican  side 
of  the  Rio  Grande.  U.S.  officials  who  pursued  three  fleeing  SUVs  to  the  Mexican  border 
saw  what  appeared  to  be  a  Mexican  military  Humvee  help  one  of  the  SUVs  when  it  got 
stuck  in  the  river. .  .When  that  didn't  work,  a  group  of  men  dressed  in  civilian  clothes 
started  unloading  what  appeared  to  be  bundles  of  marijuana  from  the  SUV,  and  the  stuck 
vehicle  was  then  torched. . .  A  second  SUV  had  a  flat  tire  and  was  left  behind  in  the 
United  States  and  its  occupant  ran  across  the  border.  (14) 

Again,  in  November  2005,  U.S.  border  patrol  agents  attempted  to  pull  over  a  suspect 
truck  on  Interstate  10  in  Texas.  The  driver,  rather  than  pulling  over,  exited  the  freeway  and  fled 
south  towards  the  Rio  Grande  and  the  Mexican  border.  While  attempting  to  cross  the  Rio 
Grande,  the  truck  got  stuck  and  the  driver  got  out  and  fled  into  Mexico  on  foot.  The  border  patrol 
found  3  tons  of  marijuana  in  the  truck  and  called  for  reinforcements  from  the  Texas  State 
Troopers.  Shortly  after  the  troopers  arrived  and  the  officers  started  emptying  the  truck,  the 
driver,  who  had  fled  into  Mexico,  returned  with  an  armed  militia  and  a  bulldozer  to  pull  the  truck 
out.  The  U.S.  Border  Patrol  and  the  Texas  State  Troopers,  outnumbered  by  the  heavily  armed 
Mexican  militia,  were  forced  to  allow  the  Mexicans  to  retrieve  the  truck,  still  two  thirds  full  of 
marijuana,  and  take  it  back  into  Mexico.  (15) 

In  yet  another  incident  “Chief  Deputy  Mike  Doyal  of  the  Hudspeth  County  [Texas] 
Sheriffs  Department  said  that  Mexican  army  personnel  had  several  mounted  machine  guns  on 
the  ground  more  than  200  yards  inside  the  U.S.  border”  (14).  Even  more  terrifying  is  the  recent 
Department  of  Homeland  Security  report  stating  that  Mexican  troops  (or  armed  paramilitary 
forces  dressed  in  Mexican  army  uniforms)  have  entered  the  U.S.  at  least  216  times  in  the  past  9 
years.  (16) 
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Because  of  the  serious  nature  of  the  situation  on  its  southern  border,  it  has  become  clear 


that  the  U.S.  must  do  more  to  protect  its  borders  and  its  citizens.  In  addition  to  hiring  more 
border  agents  (or  putting  the  National  Guard  on  the  border)  and  providing  them  with  proper 
training  and  equipment,  the  US  has  the  ability  to  help  its  border  agents  by  using  a  number  of 
high-tech  devices  to  detect  and  track  illegal  crossings  into  our  country.  Smart  fences,  multiple 
types  of  cameras,  radar  towers,  and  seismic  sensors  linked  through  a  system-wide  wireless 
communication  network  can  be  used  to  detect  and  track  subjects.  Making  use  of  these  new 
technologies,  border  agents  will  be  able  to  perform  their  jobs  better,  more  efficiently,  and  with 
higher  success  rates  than  ever  before. 

Problem  Statement 

The  United  States  has  approximately  7000  miles  of  border  with  Canada  and  Mexico.  Of 
these  7000  miles  of  border,  most  of  the  Canadian,  and  large  parts  of  the  Mexican  border,  are 
almost  completely  unprotected. 

The  U.S.  Border  Patrol  has  20  sectors  responsible  for  detecting,  interdicting  and 
apprehending  those  who  attempt  to  illegally  enter  or  smuggle  people,  including  terrorists, 
or  contraband,  including  weapons  of  mass  destruction,  across  U.S.  borders  between 
official  ports  of  entry.  (17) 

For  example,  the  El  Centro  Sector  covers  the  Riverside  and  Imperial  counties  in  southern 
California.  (17) 

Although  the  Border  Patrol  does  have  agents  assigned  to  each  section  of  the  border,  the 
available  manpower,  given  the  magnitude  of  the  task,  is  insufficient  to  cover  every  possible  entry 
point.  In  addition,  up  until  the  present  day,  the  use  of  technology  on  the  border  to  help  with 
interdiction  efforts  has  been  limited.  While  technology  is  being  used  on  some  sections  of  the 
border  presently,  what  is  needed  is  a  systematic  effort  to  implement  technological  solutions  into 
the  areas  of  the  border  patrol  effort  where  they  are  most  effective.  By  introducing  technological 
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solutions  on  the  US  border,  especially  in  the  areas  of  detection,  tracking,  and  communications,  it 
is  possible  to  help  agents  do  their  jobs  much  more  efficiently  and  effectively.  It  is  also  much 
cheaper  to  add  a  technological  infrastructure  to  help  agents  do  their  jobs,  than  it  is  to  hire  the 
large  number  of  additional  agents  that  will  otherwise  be  required.  (18) 

Research  Objective 

The  primary  objective  of  this  research  is  the  development  of  a  model  optimizing  the 
placement  of  electronic  sensors  on  a  border  network  given  a  pre-determined  budgetary 
constraint.  The  model  is  capable  of  handling  multiple  sensor  types  which  are  placed  together  as 
packages.  Also,  some  sensors  operate  during  daytime,  others  operate  during  nighttime,  and  yet 
others  operate  both  during  daytime  and  nighttime.  The  probability  of  an  intruder  being  detected, 
by  each  sensor  type,  is  calculated  for  each  node  in  the  network.  Then,  the  probabilities  for  each 
sensor  type  are  combined  using  the  assumption  of  independent  probabilities.  A  separate 
probability  of  detection  at  each  node  is  calculated  for  daytime  and  nighttime  and  both  (by  taking 
the  average  of  the  daytime  and  nighttime  probabilities).  The  model  then  uses  several  techniques 
to  place  sensors  at  nodes  in  order  to  maximize  coverage  (probability  of  detection)  on  the 
network. 

Research  Focus 

The  research  is  focused  on  creating  a  proof-of-concept  model  for  placement  of  a 
distributed  network  of  unattended  electronic  sensors  in  order  to  maximize  the  probability  of 
detecting  intruders.  The  model  maximizes  the  probability  of  detecting  illegal  aliens  using 
heuristic  methods  to  place  electronic  sensors  creating  an  interior  surveillance  network  capable  of 
detecting  intruders  after  they  have  already  breached  the  border.  The  model  will  not  account  for 
technology  placed  on  the  border  itself  (such  as  smart  fences).  As  an  additional  requirement,  the 
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model  will  be  easy  to  implement  and  modify  giving  the  user  the  ability  to  quickly  make  changes 
and  re-run  the  model  in  order  to  adapt  to  changing  requirements.  Microsoft  Excel  is  the  software 
of  choice  for  this  research  because  of  its  high  world-wide  market  share.  In  fact,  Microsoft  “owns 
more  than  90  percent  of  the  office  productivity  application  market”  through  their  Office  software 
suite,  which  contains  Excel  (19).  Therefore,  it  is  not  unreasonable  to  assume  that,  in  any 
organization,  there  is  at  least  one  person  who  knows  how  to  use  Excel;  making  the  model  much 
more  likely  to  be  used. 

Overview 

The  remainder  of  the  document  has  a  review  of  prior  interdiction  related  research  as  well 
as  the  software  which  will  be  required  to  complete  this  research.  This  is  presented  in  Chapter 
2.Then,  the  model,  which  is  object  of  this  research,  is  developed  and  tested.  Finally,  future 
research  recommendations  are  made. 
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II.  Literature  Review 


National  Border  Patrol  Strategy 

The  Border  Patrol’s  strategy  for  protecting  the  national  borders,  as  stated  in  the  National 
Border  Patrol  Strategy  of  2004,  consists  of  the  following  five  objectives: 

1 .  Establish  substantial  probability  of  apprehending  terrorists  and  their  weapons  as  they 
attempt  to  enter  illegally  between  the  ports  of  entry.  (13:7-1 1) 

2.  Deter  illegal  entries  through  improved  enforcement.  (13:7-1 1) 

3.  Detect,  apprehend,  and  deter  smugglers  of  humans,  drugs,  and  other  contraband. 
(13:7-11) 

4.  Leverage  “Smart  Border”  technology  to  multiply  the  effect  of  enforcement 
personnel.  (13:7-1 1) 

5.  Reduce  crime  in  border  communities  and  consequently  improve  quality  of  life  and 
economic  vitality  of  targeted  areas.  (13:7-1 1) 

The  Border  Patrol  has  identified  four  approaches  they  will  use  to  achieve  the  outlined 

objectives: 

1.  A  more  flexible,  well-trained,  nationally-directed  Border  Patrol.  (13:7-11) 

2.  Specialized  teams  and  rapid-response  capabilities.  (13:7-11) 

3 .  Intelligence-driven  operations .  (13:7-11) 

4.  Infrastructure,  facility,  and  technology  support.  (13:7-11) 

Ninety  percent  of  arrests  made  by  the  Border  Patrol  each  year  occur  along  the  2000  mile 
long  U.S.  border  with  Mexico. 

The  Border  Patrol  has  experienced  success  in  gaining  operational  control  of  the 
[Southern]  border  in  some  of  the  highest  trafficked  areas,  such  as  San  Diego  [CA],  El 
Paso  [TX],  and  McAllen  [TX],  However,  many  other  areas  along  the  southwest  border 
are  not  yet  under  operational  control,  and  the  daily  attempts  to  cross  the  border  by 
thousands  of  illegal  aliens  from  countries  around  the  globe  continue  to  present  a  threat  to 
U.S.  national  security.  (13:5) 
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The  Border  Patrol  has  identified  the  following  strategies  for  controlling  the  U.S. -Mexico 
(Southern)  border: 

1 .  Deter  or  deny  access  to  urban  areas,  infrastructure,  transportation,  and  routes  of 
egress  to  smuggling  organizations  through  checkpoints,  intelligence-driven  special 
operations,  and  targeted  patrols;  (13:15-16) 

2.  Expand  control  through  increased  and  more  mobile  personnel  and  improved  air  and 
ground  support;  (13:15-16) 

3.  Increase  rapid  response  capabilities;  (13:15-16) 

4.  Continue  and  expand  the  appropriate  mix  of  improved  infrastructure  and  technology; 
(13:15-16) 

a.  Sensing  systems,  Remote  Video  Surveillance  and  Sensing  (RVSS)  cameras,  air 
support,  and  Unmanned  Aerial  Vehicles  (UAVs)  (13:15-16) 

b .  Radiation  detection  equipment  (13:15-16) 

c.  Improved  communication  infrastructure  (Land  Mobile  Radio,  cellular  coverage, 
satellite  communication  capabilities)  (13:15-16) 

d.  Remote  access  to  national  law  enforcement  databases  through  the  use  of  mobile 
computing  solutions  (13:15-16) 

Network  Interdiction 

Network  Interdiction  involves  a  network  user  trying  to  utilize  a  network  to  optimize  the 
movement  of  goods  and  information,  while  a  network  interdictor  attempts  to  stop  or  reduce  the 
movement  of  material  and  information  through  the  network.  From  a  military  perspective, 
interdicting  is  generally  modeled  by  destroying  the  nodes  of  a  network  or  reducing  their 
effectiveness  below  a  predetermined  threshold.  However,  for  border  interdiction,  the  goal  is  to 
optimize  coverage  over  a  given  network  in  order  to  improve  the  success  rate  of  interdiction 
efforts.  Destroying  the  nodes  of  the  network  is  generally  not  an  option. 


Shortest  Path  Network  Interdiction. 


In  their  article,  Eitan  and  Wood  develop  a  new  method  for  maximizing  the  shortest-path 
between  two  nodes  in  a  network.  If  the  interdictor  had  an  unlimited  budget,  he  would  simply 
solve  the  cut-set  problem  and  destroy  all  of  the  designated  arcs  thus  making  it  completely 
impossible  for  the  user  to  move  anything  across  the  network.  Of  course,  this  is  rarely  the  case.  In 
reality,  there  will  be  budgetary  constraints  which  the  interdictor  must  follow.  Therefore,  while 
the  interdictor  may  not  be  able  to  completely  cut  the  network,  he  can  maximize  the  length  of  the 
shortest  path.  The  object  is  to  destroy  (or  lengthen)  a  select  number  of  arcs  in  order  to  optimize 
the  disruption  to  the  network  under  the  budgetary  constraints  placed  on  the  interdictor.  (20) 

The  shortest-path  network  interdiction  problem  can  be  solved  using  a  branch  and  bound 
plus  linear  programming  relaxation  approach.  However,  this  method  can  be  very  time 
consuming,  especially  when  dealing  with  large  networks.  Additionally,  in  the  military  realm,  the 
need  for  solutions  is  often  time  sensitive.  For  this  reason,  and  others,  the  authors  have  developed 
an  algorithm  that  improves  on  the  efficiency  of  a  linear  relaxation  solution.  (20) 

Eitan  and  Wood  started  by  formulating  the  “Maximizing  the  Shortest  Path”  (MXSP) 
problem  as  a  Mixed  Integer  Programming  (MIP)  problem.  They  also  developed  four  separate 
decomposition  methods  which  solve  problems  quicker  than  the  traditional  branch-and-bound 
linear  programming  approach.  On  its  own,  Benders’  decomposition  performed  quite  poorly,  but 
with  the  addition  of  “supervalid  inequalities”,  it  showed  significant  improvements  in 
computational  efficiency.  (20:97) 

The  “Supervalid  Inequality”  (SVI)  introduced  in  the  article  can  be  viewed  as  a 
generalized  version  of  the  “standard  valid  inequality”  (or  “cut”).  However,  whereas  the  standard 
valid  inequality  would  not  reduce  the  number  of  feasible  solutions,  the  SVI  does  indeed  reduce 
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the  number  of  feasible  solutions.  Furthermore,  the  feasible  solutions  are  reduced  in  such  a  way  as 
to  guarantee  that  the  optimal  solution  is  not  removed  (unless  the  incumbent  solution  is  also  the 
optimal).  (20:100) 

Out  of  the  four  decomposition  algorithms  developed  in  the  article,  two  work  quite  well 
with  the  MXSP  problem.  In  fact,  all  of  them  offer  an  increase  in  efficiency  over  the  classic 
branch-and-bound  LP.  However,  the  most  intriguing  aspect  of  this  article  is  the  second 
decomposition  algorithm  because  it  can  be  generalized  to  other  network  and  system  interdiction 
problems.  Indeed,  the  authors  claim  this  algorithm  is  already  being  used  to  solve  a  “tri-level 
system  defense  problem”  in  order  to  “harden  a  road  network  against  attack.”  (20: 1 10) 

The  SVIs  developed  in  this  article  proved  to  be  a  very  useful  tool.  Using  SVIs,  optimality 
was  determined  significantly  faster  than  with  Benders’  decomposition.  The  main  shortcoming  of 
all  of  these  algorithms  is  reduced  flexibility  to  one  degree  or  another.  SVIs  are  an  effective  tool 
for  more  efficient  and  faster  solutions,  but  they  can  only  be  used  for  a  reduced  set  of  problems.  If 
time  is  not  of  the  essence,  it  may  be  easier  and  simpler  to  employ  the  classic  branch  and  bound 
plus  LP  relaxation  technique  which  will  theoretically  solve  all  IPs  and  MIPs  eventually.  (20) 

LP  Optimization. 

Pulat  [2005]  develops  a  mixed  integer  linear  program  which  optimizes  border  interdiction 
in  the  Yuma  sector  of  the  U.S. -Mexico  border.  He  studies  scenarios  where  the  intruder  is 
traveling  in  a  vehicle  and  scenarios  where  the  intruder  is  traveling  on  foot.  The  scenarios  are 
further  divided  into  the  case  where  the  intruder  knows  the  U.S.  Border  Patrol’s  positions  ahead 
of  time,  versus  the  case  where  the  intruder  is  not  pre-aware  of  these  positions.  Pulat  also  makes  a 
“distinction  between  actions  that  can  only  lead  to  detection  [sensors,  helicopters]  and  action  that 
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can  also  lead  to  capture  in  addition  to  detection  [road  patrols,  checkpoints,  remote  observation 
posts].”  (21:39) 

Pulat  uses  a  network  representation  of  nodes,  arcs,  and  centers  of  land  parcels  overlaying 
a  satellite  map  of  the  Yuma,  Arizona  border  area.  He  uses  open  source  information  from  the 
Border  Patrol  and  identifies  all  candidate  defensive  actions  based  on  the  location  of  checkpoints, 
road  patrols,  off-road  operations,  remote  observation  posts,  and  electronic  sensors.  He  also 
identifies  intruder  actions  and  creates  a  “Two-Sided  Mixed  Integer  Optimization  Model  to 
Minimize  Maximum  Probability  of  Escape”  (21:25).  Using  a  number  of  different  scenarios,  Pulat 
identifies  critical  road  segments  and  land  parcels  to  be  defended  and  studies  the  “effects  of 
employing  different  types  of  assets  and  strategies  on  the  infiltration  patterns.”  (21:39) 

Continuous  Network  Interdiction. 

Washburn  [2006]  develops  a  network  interdiction  model  for  economic  networks  with 
indefinite  time  outlooks.  This  model  seeks  to  minimize  the  fraction  of  product  that  makes  it  from 
its  origination  point  to  its  destination  point  without  being  interdicted.  The  model  is  developed  as 
a  two-person  zero-sum  game.  He  also  explores  the  consequences  of  allowing  the  interdictor  to 
sell  confiscated  goods.  This  not  only  increases  the  interdictor’ s  budget,  leading  to  a  larger 
interdiction  effort,  but  also  depreciates  the  commodity  making  it  harder  for  the  shipper  to  make  a 
profit.  “This  leads  to  a  Nash  equilibrium  where  the  shipper’s  quantity  shipped  is  in  equilibrium 
with  the  interdictor’s  budget  for  interdiction.”  (22:1) 

Game  Theory  Approach. 

Washburn  and  Wood  [1995]  develop  a  game  theory  approach  to  network  interdiction. 

The  game  takes  place  on  a  network  of  nodes  and  arcs  with  one  evader  and  one  interdictor.  For 
each  arc  in  the  network,  a  constant  probability  of  detection  is  determined  beforehand.  Then, 


11 


while  the  evader  determines  a  “path-selection”  strategy  minimizing  the  probability  of  detection, 
the  inspector  determines  an  “arc-inspection”  strategy  maximizing  the  probability  of  detection. 
The  authors  show  that  this  type  of  problem  can  be  solved  using  standard  network  flow 
techniques.  They  also  discuss  problems  with  “unknown  origins  and  destinations”  as  well  as 
“multiple  interdictors  and  evaders.”  (23:243) 

Sensor  Placement 

Remote  sensing  technologies  have  the  potential  of  greatly  reducing  the  number  of 
personnel  needed  for  border  patrol  while  at  the  same  time  increasing  the  probability  of  detecting 
and  capturing  intruders.  While  the  border  patrol  has  been  using  a  limited  number  of  electronic 
sensors  and  other  devices  for  a  number  of  years,  they  do  not  have  an  integrated  electronic 
network  of  sensors  designed  to  detect,  track,  and  aid  in  the  capture  of  illegal  aliens  and 
smugglers. 

Sensor  Placement  Algorithms  for  Effective  Coverage. 

Dhillon  and  Chakrabarty  “present  two  algorithms  for  the  efficient  placement  of  sensors  in 

a  sensor  field.”  (24:1609)  Both  algorithms  are 

aimed  at  optimizing  the  number  of  sensors  and  determining  their  placement  to  support 
distributed  sensor  networks.  The  optimization  framework  is  inherently  probabilistic  due 
to  the  uncertainty  associated  with  sensor  detections.  The  proposed  algorithms  address 
coverage  optimization  under  the  constraints  of  imprecise  detections  and  terrain 
properties.  These  algorithms  are  targeted  at  average  coverage  as  well  as  at  maximizing 
the  coverage  of  the  most  vulnerable  grid  points.  The  issue  of  preferential  coverage  of  grid 
points  (based  on  relative  measures  of  security  and  tactical  importance)  is  also  modeled. 
(24:1609) 

For  both  algorithms,  it  is  assumed  that 

the  probability  of  detection  of  a  target  by  a  sensor  varies  exponentially  with  the  distance 
between  the  target  and  the  sensor.  A  target  with  distance  d  from  a  sensor  is  detected  by 
that  sensor  with  the  probability  e‘ad.  The  parameter  a  can  be  used  to  model  the  quality  of 
a  sensor  and  the  rate  at  which  its  detection  probability  diminishes  with  distance. 

(24:1610) 
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For  every  set  of  points  i  and  j  in  the  sensor  field,  two  probability  values  are  assigned:  py, 
which  denotes  the  probability  that  a  target  at  point  j  is  detected  by  a  sensor  at  point  i,  and  py, 
which  denotes  the  probability  that  a  target  at  point  i  is  detected  by  a  sensor  at  point  j.  The 
probabilities  py  and  p,,  are  symmetric  in  most  cases  but  can  differ  in  the  presence  of  obstacles. 

Dhillon  and  Chakrabarty’s  first  algorithm  (MAX_AVG_COV)  attempts  to  maximize  the 
average  coverage  of  the  grid  points,  while  their  second  algorithm  (MAX_MIN_COV)  attempts 
to  maximize  the  coverage  of  the  grid  point  which  is  least  effectively  covered;  that  is  the  grid 
where,  if  located,  a  target  would  have  the  least  probability  of  being  detected.  Dhillon  and 
Chakrabarty  test  the  two  algorithms,  on  an  8  by  8  grid,  against  each  other  as  well  as  random  an 
uniform  placement  of  sensors.  They  conclude  that  the  MAX_MIN_COV  algorithm  produces 
superior  results,  i.e.  they  achieved  the  best  probability  of  detection  (coverage)  using  this 
algorithm.  Furthermore,  they  discuss  continued  research  which  would  include  minimum  and 
maximum  ranges  for  each  sensor.  (24:1614) 

Sensor  Placement  Algorithm  for  Minimalistic  Grid  Coverage. 

Dhillon,  Chakrabarty,  and  Iyengar  present 

a  resource  bounded  optimization  framework  for  sensor  resource  management  under  the 
constraints  of  sufficient  grid  coverage  of  the  sensor  field.  The  proposed  theory  is  aimed  at 
optimizing  the  number  of  sensors  and  determining  their  placement. .  .The  proposed 
algorithm  addresses  optimization  under  constraints  of  imprecise  detections  and  terrain 
properties.  The  issue  of  preferential  coverage  of  grid  points  (based  on  relative  measures 
of  security  and  tactical  importance)  is  also  modeled.  (25:1) 

For  every  set  of  points  i  and  j  in  the  sensor  field,  two  probability  values  are  assigned:  py, 
which  denotes  the  probability  that  a  target  at  point  j  is  detected  by  a  sensor  at  point  i,  and  py, 
which  denotes  the  probability  that  a  target  at  point  i  is  detected  by  a  sensor  at  point  j.  The 
probabilities  py  and  py  are  symmetric  in  most  cases  but  can  differ  in  the  presence  of  obstacles. 
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The  algorithm  uses  an  iterative  “greedy  heuristic”  to  determine  the  best  placement  of  a 
single  sensor  one  at  a  time.  At  every  iteration,  the  algorithm  adds  one  sensor  and  calculates  the 
new  probabilities  for  the  entire  grid.  It  also  keeps  track  of  improvements  from  previous 
iterations.  The  algorithm  continues  placing  sensors  until  the  miss  probability  for  each  point  is 
smaller  than  the  maximum  permitted  value.  Preferential  coverage  areas  in  the  grid  can  be 
implemented  by  lowering  the  maximum  miss  probability  for  preferred  points  and  thereby  forcing 
a  higher  probability  of  detection  in  those  areas.  Also,  the  algorithm 

makes  the  implicit  assumption  that  sensor  detections  are  independent,  i.e.  if  a  sensor 

detects  a  target  at  a  grid  point  with  probability  pi,  and  another  detects  the  same  target  at  a 

grid  point  with  probability  p2,  then  the  miss  probability  for  the  target  is  (l-pl)(l-p2). 

(25:4) 

The  algorithm  presented  by  Dhillon,  Chakrabarty,  and  Iyengar  adds  one  sensor  at  a  time 
to  the  grid  until  certain  preset  conditions  are  met.  It  is  intended  to  determine  the  minimum 
number  of  sensors  needed  to  meet  the  preset  requirements.  It  does  not  backtrack  in  order  to  find 
the  optimum  placement  of  sensors  at  each  iteration. 

Sensor  Technology 

There  are  many  useful  sensor  technologies  which  can  be  employed  by  the  Border  Patrol 
for  intruder  detection.  Some  of  them  are  discussed  below. 

Cameras. 

There  are  a  large  number  of  camera  systems  and  technologies  available  from  various 
defense-focused  vendors.  These  include  the  more  traditional  daylight  cameras,  low-light  level 
cameras,  and  infra-red  (IR)  cameras,  as  well  as  the  newer  and  more  sophisticated  Forward 
Looking  Infra-Red  (FLIR)  and  Range-Gated  cameras  (26:1).  FLIR  cameras  are  thermal  imaging 
cameras.  Unlike  traditional  IR  cameras,  FLIR  cameras  do  not  require  IR  illuminators,  which 
make  it  almost  impossible  for  intruders  to  spot  them.  Unfortunately,  FLIR  cameras  do  have  some 
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significant  drawbacks  as  they  do  not  work  well  in  adverse  weather  conditions  and  they  can  be 
evaded  by  using  techniques  which  minimize  heat  signatures  (26:7).  Range-Gated  cameras  make 
use  of  lasers  and  other  technologies  to  literally  see  through  snow,  rain  and  fog  at  any  time  of  the 
day  or  night  (26:8).  Two  examples  of  camera  systems  are  included  below,  but  many  others  are 
available.  Also  note  that,  as  technology  continues  to  improve,  the  included  examples  will  be 
outdated. 

The  GVS1000  (see  Figure  1)  is  a  “long-range  active-infrared  day/night  surveillance 
system.”  It  delivers  1.2  kilometers  of  “classification  level”  surveillance  in  complete  darkness.  It 
also  contains  integrated  software  which  can  “classify,  recognize,  and/or  identify”  targets.  (1) 


Figure  1.  GVS1000  Long-Range  Surveillance  System  (1) 


The  Axsys  ExtremeXS  thermal  imaging  camera  (see  Figure  2)  is  a  “rugged  camera  with 
extensive  detection  capabilities.”  The  ExtremeXS  can  detect  a  human  sized  target  at  up  to  4.5 
kilometers  distance  in  less  than  ideal  conditions.  (2) 


15 


Figure  2.  Axsys  ExtremeXS  Thermal  Imagery  Camera  (2) 


Ground  Surveillance  Radar. 

The  Motorola  Modular  Surveillance  Radar  is  a  man  portable  radar  system  capable  of 
detecting  a  single  person  sized  object  up  to  3  miles  away.  It  can  detect  a  small  vehicle  up 
to  7  miles  away,  and  a  larger  vehicle  up  to  12  miles  away.  Groups  of  vehicles  or  people 
improve  detectability.  The  MSR  provides  target  location  accuracy  of  15  meters  in  range, 
and  .6  degrees  in  azimuth.  The  radar  system  has  been  mounted  to  vehicles,  trailers,  and 
fixed  site  towers,  and  has  been  used  operationally  since  1990.  The  radar  system  provides 
wide  area  surveillance.  When  a  daylight/infrared  camera  system  is  used  in  combination 
with  the  radar,  target  identification  is  possible.  The  radar  can  be  remotely  controlled  by 
radio  link,  or  long  haul  RS-232  lines.  (3) 

Additionally,  Dragoon  Technologies  has 

developed  a  modern  map  based  application  for  detected  target  display.  The  application 
can  be  used  to  steer  additional  sensors  and  accepts  GPS  input  for  mobile  applications. 
The  radar  system  utilizes  mil-spec  construction  and  operates  on  24VDC.”  (3) 

Figure  3  shows  the  Motorola  MSR-20  Ground  Search  Radar  mounted  on  a  tower  along 

with  a  video  camera  and  infra-red  sensor. 
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Figure  3.  Motorola  MSR-20  Ground  Search  Radar  (3) 


Seismic  Sensors. 

There  are  a  number  of  promising  seismic  sensing  technologies  which  can  be  used  for 
border  security.  Maier  [2004]  developed  a  seismic  intrusion  sensor  technology  which  uses  buried 
fiber  optic  cables,  lasers,  and  piezoelectric  transducers,  to  detect  and  locate  walking  intruders  at 
distances  up  to  2  kilometers  away  (37).  In  2002,  the  U.S.  Army  Corps  of  Engineers  tested  a 
different  approach  to  a  seismic  sensing  system.  This  system  used  nodes  made  up  of  six  sensors 
placed  in  a  circular  area  with  a  6  meter  diameter.  The  test  concluded  that  seismic  sensors  were 
effective  at  distances  up  to  1  km  away  under  simulated  battlefield  environments.  (38) 

Palm  PDA  Based  Intelligence  Distribution. 

Getting  sensor  data  collected,  processed,  and  distributed  to  officers  in  the  field  can  be  a 

lengthy  process  if  it  involves  human-in-the-loop  interactions.  Dragoon  Technologies  has 

developed  two  applications  [Figure  2]  for  PDA  computers  that  put  MTI  data,  video  and 
freeze  frame  imagery  onto  a  computer  slightly  larger  than  a  deck  of  playing  cards.  The 
screens  are  sunlight  viewable  and  the  form  factor  is  soldier/operative  friendly. 
Communications  to  these  devices  is  currently  available  as  RS-232,  TCP/IP,  WiFi, 
Bluetooth,  and  cellular  telephone.  The  links  all  include  the  ability  to  send  bi-directional 
data  to  include  GPS  position  of  the  Palm  PDA  to  a  server  and  chat  messaging.  Five 
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streaming  video  is  now  available  as  well.  The  PDA  represents  the  future  of  intelligence 
distribution  to  scaled,  mobile  devices.  (4) 


NOT  data  as  nest  ed  on  Palin  PDA 


Imagery  zoomified  and  viewed  on  Palm  PDA 


Figure  4.  Palm  PDA  based  intelligence  distribution  (4) 


PDA  devices  have  the  potential  to  provide  border  agents  in  the  field  with  near 
instantaneous  information  enabling  them  to  track  and  capture  intruders  with  unprecedented 
ease. 


Software 

In  addition  to  Microsoft  Excel,  a  few  other  software  programs  are  needed  to  create  the 
proof-of-concept  network.  Since  the  Border  Patrol  does  not  make  their  maps  and  mapping 
software  available  to  the  public,  Google  Earth  is  used  as  the  base  map  for  the  network. 

Google  Earth. 

Google  Earth  is  a  free  virtual  globe  mapping  program  originally  developed  by  Keyhole 
Inc.,  but  later  purchased  and  distributed  by  Google  Inc.  It  provides  satellite  images  of  the  entire 
world  with  overlaid  road  maps  and,  in  some  places,  3D  terrain  and  building  models.  The 
program  also  allows  users  to  create  and  store  their  own  points  of  interest  called  “place  marks”. 
(6) 
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GEPath. 


GEPath  is  a  freeware  program  “developed  to  make  paths  and/or  draw  circles  and 
polygons  with  place  marks  saved  by  Google  Earth.”  It  parses  Google  Earth  “kml”  files  (kml  files 
are  files  written  in  extensible  markup  language  and  used  by  Google  Earth  to  show  user  specified 
information)  and  retrieves  place  mark  information  such  as  the  place  mark’s  name,  latitude,  and 
longitude. 

The  data  can  also  be  typed  into  the  application  or  pasted/exported  to  the  clipboard.  Files 
generated  by  GE-Path  are  exported  to  Google  Earth.  This  application  calculates 
distances,  bearing  and  area.  (7) 

Frontline  Systems  OptQuest  Solver. 

The  Excel  Solver  allows  the  user  to  optimize  a  given  objective  function  based  on  a  set  of 
changeable  cells  (variables)  and  a  set  of  constraints.  Frontline  Systems  is  the  company  which 
developed  the  Excel  Solver  for  Microsoft.  However,  the  Excel  Solver  is  limited  in  its  scope.  It 
has  a  maximum  capacity  of  200  variables  and  constraints  for  linear  models  and  can  only  solve  a 
limited  number  of  non-linear  models.  Frontline  Systems  offers  the  Premium  Solver  Platform  and 
a  number  of  “field  installable  engines”  to  extend  the  capability  of  the  Excel  Solver.  Specifically, 
the  OptQuest  solver  (one  of  the  field-installable  engines)  “employs  metaheuristics  such  as  tabu 
search  and  scatter  search  to  solve  nonsmooth  optimization  problems  of  up  to  5,000  variables  and 
1,000  constraints.  It  also  supports  integer  variables.”  While  not  guaranteeing  an  optimal  solution, 
the  OptQuest  Solver  “finds  remarkably  good  solutions  with  unprecedented  speed.”  (27) 

In  the  next  chapter.  Border  Patrol  interdiction  and  the  need  for  innovations  is  discussed  in 
detail.  After  explaining  the  need  for  technological  innovation  for  border  security,  a  model  is 
developed  optimizing  the  placement  of  electronic  sensors  in  order  to  maximize  the  probability  of 
detecting  intruders. 
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III.  Methodology 


In  this  chapter,  the  traditional  approach  to  border  patrol  is  discussed  along  with  the  need 
for  a  new  approach.  Then,  a  model  is  built  for  placement  of  distributed  sensors  on  a  network 
with  the  goal  of  maximizing  the  probability  of  detecting  intruders.  This  model  is  intended  to  be 
the  first  part  of  the  overall  strategy  of  creating  a  new  technological  approach  to  border  patrol. 

Traditional  approach 

Traditionally,  border  protection  has  been  a  very  manpower  intensive  job.  The  job  requires 
many  border  patrol  agents  in  vehicles,  on  horseback,  or  on  foot  to  patrol  areas  searching  for 
intruders.  Intruder  detection  can  also  be  performed  by  helicopter  patrols,  but,  while  helicopters 
greatly  improve  speed  and  the  probability  of  detection,  they  are  expensive  to  purchase,  fly,  and 
maintain.  Once  intruders  are  detected,  the  patrol  agents  must  change  tasks  and  attempt  to 
apprehend  the  intruders. 

Another  traditional  method  for  border  protection  is  the  interior  checkpoint.  The  Border 
Patrol  uses  both  permanent  and  temporary  immigration  checkpoints  where  all  vehicle  traffic  is 
stopped  in  order  to  detect  and  apprehend  illegal  aliens,  drugs,  and  other  illegal  activity.  The 
permanent  checkpoints  are  generally  located  on  national  roads  and  interstate  freeways,  while 
temporary  checkpoints,  called  “tactical  checkpoints,”  are  located  on  smaller  arterial  and  rural 
streets  with  traffic  volumes  as  small  as  a  few  hundred  vehicles  per  day  (5).  The  2005 
Government  Accountability  Office  report  on  immigration  checkpoints  claims  that  “while 
changing  locations  of  tactical  checkpoints  would  appear  to  offer  the  potential  element  of 
surprise. . .  the  border  patrol  [claims]  that  the  smugglers  of  aliens  and  contraband. .  .use  cell 
phones  and  communications  equipment  to  alert  confederates  of  the  presence  of  checkpoints 
within  minutes  of  their  being  relocated”  (5:23,24).  However,  despite  the  fact  that  smugglers  have 
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become  increasingly  more  sophisticated  in  their  use  of  technology,  there  is  sufficient  reason  to 
believe  that  checkpoints  make  up  a  useful  part  of  a  multi-layered  border  protection  strategy.  For 
example,  in  2004,  the  Border  Patrol’s  Southwest  interior  checkpoints  used  10  percent  of  the 
region’s  border  patrol  agents,  contributed  to  8  percent  of  the  total  number  of  apprehensions,  and 
31  percent  of  marijuana  and  74  percent  of  all  cocaine  seizures.  (5:29,30) 


Figure  5.  Tactical  Immigration  Checkpoint  (5:29) 


Checkpoints  are  generally  effective  only  against  vehicular  traffic  because  pedestrians 
tend  to  find  ways  around  them.  However,  if  strategically  placed,  it  is  possible  for  checkpoints  to 
act  as  temporary  deterrents  against  pedestrian  intruders.  (5) 

The  Need  for  Technological  Innovation 

With  half  a  million  people  crossing  the  border  illegally  each  year,  it  is  obvious  that  the 
border  patrol  is  not  able  to  stop  all  of  the  illegal  cross-border  inflow  of  aliens,  drugs,  and  other 
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contraband.  It  is  felt  that  the  U.S.  Border  Patrol  is  undermanned  and  underfunded  but 

Washington  has  done  little  to  change  this  situation;  even  after  the  events  of  September  11,  2001 

(28).  In  May  of  2006,  President  Bush  announced  a  $1.9  billion  plan  which  has  placed  nearly 

6000  National  Guard  troops  to  the  U.S.  border  with  Mexico  (29).  Unfortunately,  National  Guard 

troops  now  stationed  on  the  Mexican  border  cannot  be  fully  utilized  because 

under  existing  rules  of  force  signed  by  the  Department  of  Defense  and  border  state 
governors,  soldiers  are  not  supposed  to  stop,  arrest,  or  shoot  armed  illegal  immigrants. 
They  are  instructed  only  to  look,  listen  and  report  their  location  to  the  Border  Patrol.  (30) 

While  putting  the  National  Guard  on  the  border  may  be  a  great  idea,  ordering  the  Guard  to 

maintain  the  status  of  observers  turns  them  into  nothing  more  than  a  human  sensor  network.  This 

job  can  be  done  more  effectively,  and  possibly  cheaper,  with  an  electronic  sensor  network. 

As  with  almost  all  organizations,  the  largest  part  of  the  Border  Patrol’s  budget  goes  to 

payroll.  This  makes  it  very  difficult  to  add  additional  manpower  because  it  requires  a  large 

budgetary  increase.  In  fact,  even  if  the  Border  Patrol  was  appropriated  enough  funds  to  double  its 

manpower,  it  would  not  guarantee  significantly  better  results.  After  all,  the  Border  Patrol  is 

currently  only  able  to  capture  an  estimated  25%  of  intruders  (12).  Even  if  the  Border  Patrol 

managed  to  cut  down  the  rate  of  illegal  border  crossings  to  half,  or  even  one-fourth,  of  their 

current  rate,  there  would  still  be  a  serious  illegal  immigration  problem. 

This  is  where  technological  innovation  can  be  used  as  a  force  enhancer.  Installing  smart 

fences  on  the  nation’s  borders  would  allow  agents  to  know  exactly  when  and  where  a  breach 

occurs.  Installing  sensor  packages,  including  radar,  video  cameras,  infrared  cameras,  seismic 

sensors,  and  other  advanced  technologies,  would  allow  agents  to  detect  and  track  intruders  in  real 

time.  This  would  eliminate  the  most  time  consuming  part  of  an  agent’s  job  (searching  for 


22 


intruders)  and  allow  the  agent  to  focus  most,  if  not  all,  of  his  or  her  time  on  apprehending 
intruders. 

Thus,  rather  than  just  adding  more  agents,  it  is  essential  that  the  Border  Patrol  provide  its 
agents  with  the  latest  advanced  technology  to  help  them  do  their  jobs  safely  and  much  more 
effectively.  In  fact,  properly-employed  technology  acts  as  a  force  multiplier  for  border  security 
personnel.  (13) 

New  Approach 

While  purchasing  new  sensors  and  other  technologies  for  the  border  patrol  is  very 
important,  the  funds  will  be  less  effective  if  the  new  technology  is  not  properly  employed.  Given 
a  set  budget,  it  is  extremely  important  that  the  border  patrol  be  able  to  balance  the  training  and 
sustainment  of  personnel  as  well  as  technology  and  infrastructure  (13).  The  Border  Patrol  must 
be  able  to  identify  how  many  sensors  they  need  to  buy  and  where  they  should  place  them,  based 
on  reasonable  budgetary  constraints.  A  computer  software-focused  approach  will  be  employed  to 
help  the  Border  Patrol  make  this  important  decision. 

It  was  decided  to  use  a  network  manually  created  and  overlaid  onto  a  map  using  Google 
Inc.’s  free  Google  Earth  software.  The  node  coordinates  are  imported  into  Excel  and  used  to 
populate  an  optimization  model.  The  model  is  created  using  various  techniques  to  optimize  the 
purchase  and  placement  of  electronic  sensors  under  pre-determined  budgetary  constraints. 

Data  Development 

The  development  of  the  network  was  done  in  several  steps.  First  a  location  was  picked 
for  the  network.  Then,  the  nodes  of  the  network  were  overlaid  on  a  map  of  the  network  location. 
Finally,  the  coordinate  locations  of  the  nodes  were  extracted  and  imported  into  Excel. 
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Location. 


A  20  kilometer  section  of  the  U.S. -Mexico  border  near  Calexico,  California  was  chosen 
as  the  location  for  the  sensor  network.  Overhead  satellite  imagery  provided  by  Google  Earth 
suggests  the  area  is  comprised  almost  entirely  of  level  farmland,  with  a  uniform  elevation  and 
few  obstacles.  However,  aside  from  the  overhead  satellite  image,  little  else  is  known  about  the 
location.  In  lieu  of  a  thorough  on-ground  inspection  of  the  location,  the  network  created  from  the 
image  is  treated  as  a  notional  network.  The  assumptions  that  have  been  made  about  this  network 
may  or  may  not  represent  the  actual  conditions  at  the  location. 

Google  Earth. 

The  database  for  the  network  was  created  using  Google  Earth’s  “placemark”  feature  with 
a  simple  circle  and  diamond  node  representation. 


Figure  6.  Intersection  nodes  and  centers  of  land  parcels  (6) 

Circles  were  used  to  indicate  an  intersection  between  two  or  more  roads,  while  diamonds 
were  used  to  represent  centers  of  land  parcels.  Because  the  land  near  Calexico  is  mostly 
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agricultural,  there  was  a  need  to  differentiate  between  the  two  types  of  nodes  (i.e.  the  nodes 
indicating  centers  of  land  masses  is  only  used  for  detection  and  capture  of  walking  intruders). 

Data  Conversion. 

Google  Earth  saves  user-created  data  points  in  a  xml  document  called  a  kml  file. 
Although  kml  files  are  text  files,  their  format  makes  it  difficult  to  directly  import  their  contents 
into  Excel.  The  problem  arises  from  the  fact  that  kml  files  are  written  in  extensible  markup 
language  (xml)  and  contain  a  number  of  rows  for  each  node  in  the  network.  These  rows  contain 
xml  tag  information  as  well  as  the  node  coordinates.  The  software  program  GEPath  was  used  to 
parse  the  Google  Earth  kml  file  and  create  a  simple  spreadsheet  grid  with  the  number  of  each 
node,  its  longitude,  and  its  latitude.  (7) 


S3  GEPath  (v  1.4.1) 


File  Run  Help 

Place  Names 


P  [Add  sequential  numbers  to  names; 
I-  Append  coordinates  to  names 
P  Remove  names  and  icons 
Altitude  mode:  Clamped  to  ground  ▼ 


Make  Path  (  Loc 

p  Make  path  I-  Link  last  to  first  ^  ^pp 

Line  color  P  Arrowed  Make 

Line  width  77"/ .  E~°l  T  Mai 

Altitude  mode:  Clamped  to  ground  d  Length 


Place  Descriptions 


(•  Don't  change 
C  None 
C  Coordinates 
C  Distance  to  next  place 
C  Route  to  next  place 
C  Distance  and  route 
C  Coordinates,  distance  and  route 
P  Append  to  existing  description 


Make  Polygons 

P  Make  polygons  Segments 
Length  |  j  20  jr| 


Make  Grid 

I-  Make  grid  Line  color 
Distance  between  lines  Line  width  |  10  | 

Number  of  lines  Altitude  mode:  Clamped  to  ground  ▼ 

Area/Perimeter  (of  the  path's  projection  in  a  flat  earth  surface) 

P  Calculate  perimeter  |  |  Lenght  unit:  Kilometer 

|“  Calculate  area  |  |  Area  unit:  Kilometei2 

Decimal  places  2  P  Create  Placemark  for  area 

P  Don't  save  the  kml  file.  Show  area  and  perimeter  only  p  Open  in  GE 


P  Outline 
P  Fill 


|  Lenght  unit:  Kilometer 

|  Area  unit:  Kilometei2 

P  Create  Placemark  for  area 


801 

32.746129 

-115.594019 

0 

800 

32.617501 

-115.582621 

0 

673 

32.672540 

-115.674096 

0 

672 

32.674581 

-115.689786 

0 

671 

32.677046 

-115.690149 

0 

670 

32.675987 

-115.685055 

0 

669 

32.675233 

-115.678384 

0 

668 

32.676637 

-115.678385 

0 

667 

32.680848 

-115.677937 

0 

666 

32.684300 

-115.678346 

0 

665 

32.688485 

-115.678809 

0 

664 

32.692116 

-115.679476 

0 

663 

32.695468 

-115.679810 

0 

662 

32.699219 

-115.679102 

0 

661 

32.704885 

-115.679297 

0 

660 

32.709713 

-115.679917 

0 

659 

32.713221 

-115.686086 

0 

658 

32.713185 

-115.680679 

0 

657 

32.712828 

-115.676268 

0 

656 

32.711674 

-115.671572 

0 

ecG  rid,  www.  encoreconsultinq.  com.  au 
Clear  grid  I  Paste  grid  from  clipboard 


Delete  selected  row  | 

N. places:  674 

Repeat  rowl  altitude  | 

Figure  7.  GEPath  with  data  from  the  Calexico  kml  file  (7) 
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Microsoft  Excel  Import. 

The  grid  created  by  GEPath  (Figure  7),  was  copied  into  an  Excel  document.  Once  copied 
into  Excel,  it  became  the  foundation  for  the  sensor  placement  model. 
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Figure  8.  Data  imported  into  Excel 


Model  Development 

Using  the  data  imported  from  Google  Earth,  an  iterative  sensor  placement  model  was 
developed.  The  model  maximizes  the  probability  of  detecting  intruders  by  optimizing  the  build¬ 
up  of  a  distributed  sensor  network  subject  to  a  budgetary  constraint.  Several  different 
optimization  algorithms  are  developed  for  use  with  the  model.  Additionally,  the  model  is 
compatible  with  the  commercial  OptQuest  solver  software. 

Variables. 

For  each  of  the  673  nodes  in  the  Calexico  network,  there  is  a  binary  variable  i: 1-673, 
which  is  used  to  select  the  location  of  sensors  on  the  network.  st  is  equal  to  one  at  nodes  with 
sensors  and  zero  at  nodes  without  sensors. 
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Distances  between  Nodes 


A  distance  matrix,  dy,  was  created  calculating  the  distance  between  all  pairs  of  nodes 
(673x673)  in  the  network. 
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Figure  9.  Distance  matrix 

The  distances  were  calculated  using  the  Great  Circle  Distance  formula.  (3 1 ;  32) 

n(/?  *  co)  +  cos  (a  *  co)  *  cos(/?  *  co) 

*  xp 


/?  =  Longitude  of  Node  i 
v  =  Longitude  of  Node  j 
xp  =  1.8520  (kilometers  per  nautical  mile) 

Note  that,  while  the  Calexico  network  may  not  require  the  use  of  the  Great  Circle 
Distance  formula  (rectilinear  calculations  could  have  been  used  due  to  the  relatively  short 
distances  involved),  the  Great  Circle  Distance  formula  was  used  in  order  to  provide  scalability  to 


Sij  =  60  *  acos(sin(a  *  co)  * 

*  cos((/c  —  v)  *  co)) 

x  =  arc  length 

a  =  Latitude  of  Node  i 

/ 1  =  Latitude  of  Node  j 
n 
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the  model.  In  addition,  the  complete  grid  of  673x673  distances  is  calculated  only  one  time  and 
the  model  is  not  encumbered  by  distance  calculations. 

Sensor  Locations. 

The  model  assumes  sensors  are  placed  together  in  packages.  Each  node  has  a  defined 
package  which  may  contain  all  or  some  of  the  sensors.  The  binary  input  parameters  Sib  i: 1-673, 
k:l-4  describe  which  sensor  types  ( k )  can  be  placed  at  each  node.  i.e.  s-,k=  1  if  the  node  i  can  host 
sensor  type  k  and  Sik=0  if  the  node  i  cannot  host  sensor  type  k.  These  inputs  to  the  model  are 
made  based  on  geographical,  political,  and  economic  considerations  and  vary  based  on  the 
location  of  the  network.  For  the  network  tested  in  this  research,  it  was  decided  not  to  place  any 
seismic  sensors  at  intersection  nodes  because  their  effectiveness  to  detect  intruders  on  foot  will 
likely  be  compromised  by  legal  vehicular  traffic.  Once  the  packages  are  determined,  a  sensor 
selection  (si)  at  a  node  selects  all  sensor  types  available  to  that  node.  This  assumption  creates 
fewer  physical  sensor  locations  making  it  easier  to  secure  and  cheaper  to  maintain  the  network 
than  if  each  sensor  is  allowed  its  own  location.  In  addition,  by  placing  sensors  in  packages,  the 
number  of  variables  (si)  is  limited  to  673  regardless  of  the  number  of  sensor  types  being  used. 

Sensor  Ranges  and  Probability  Distributions. 

A  review  of  sensor  placement  literature  has  revealed  a  couple  of  different  methods  used 
to  define  a  sensor’s  probability  of  detection.  The  first,  an  unbounded  method,  uses  a  parameter  a 
to  obtain  a  probability  of  detection  of  an  intruder  by  a  sensor  which  varies  exponentially  with  the 
distance,  d,  between  the  intruder  and  the  sensor.  Using  this  method,  the  probability  of  detection 
becomes  e'aS  (24:1610).  The  second  method  is  bounded,  but  assumes  a  binary  probability  of 
detection  (d)  so  that  d=l  when  the  intruder  is  within  range  of  the  sensor  and  d=0  when  the 
intruder  is  out  of  range  of  the  sensor.  (25:3) 
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For  the  current  model,  a  new  bounded  method  for  determining  probability  of  detection  is 
developed.  This  method  places  lower  and  upper  bounds  on  each  sensor’s  probability  of  detection 
and  assumes  a  continuous  distribution  between  those  bounds.  The  cumulative  distribution 
function  (cdf)  of  the  Beta  probability  distribution  is  used  to  model  the  detection  probability  curve 
for  each  sensor  within  its  prescribed  range.  The  Beta  probability  distribution  was  chosen  for  its 
extreme  versatility.  By  changing  the  distribution’s  shape  parameters  ( a  and />'),  the  beta  density 
function  can  be  decreasing,  increasing,  convex,  concave,  uniform,  and  so  forth.  For  a  more 
detailed  description  of  how  the  parameters  a  and  [:>  affect  the  shape  of  the  Beta  probability 
distribution,  see  Appendix  A.  The  flexibility  provided  by  the  Beta  distribution  allows  the  user  to 
change  the  parameters,  and  therefore  the  curve  of  the  probability  distribution,  to  match  that  of 
the  sensors  being  used. 

The  Beta  distribution  is  available  as  a  function  in  Excel  and  requires  5  input  parameters: 

S  =  distance  to  be  evaluated,  a  and  /?  are  shape  parameters,  and  a  and  b  are  the  lower  and  upper 
bounds.  The  Beta  cdf  is  equal  to  0  at  the  lower  bound  and  1  at  the  upper  bound.  The  Beta  cdf 
will  be  assumed  to  indicate  the  probability  of  a  miss  (m)  with  0%  chance  of  a  miss  at  the  lower 
bound  (set  to  zero)  and  100%  chance  of  a  miss  at  the  upper  bound  or  beyond  (set  to  the  range  of 
each  sensor  type).  Note  that,  while  this  is  a  reasonable  assumption,  it  is  notional  and  has  not  been 
validated  from  actual  sensor  data.  It  is  assumed  that  sensors  have  a  100%  chance  of  detection  at  a 
distance  of  0  kilometers,  and  a  0%  chance  of  detection  at  a  given  distance,  with  decreasing 
probability  of  detection  between  the  given  bounds.  In  order  to  obtain  a  probability  of  detection 
(d)  equal  to  1  at  the  lower  bound  (shortest  distance)  and  0  at  the  upper  bound  (longest  distance), 
d=l-  m  is  used. 
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Probability  of  Missing  (1-  Probability  of  Detection). 

For  every  set  of  nodes  i  and  j  in  the  sensor  field,  and  for  every  sensor  type  k,  the 
probability  m*,*,  which  denotes  the  probability  that  a  target  at  node  j  is  missed  by  a  sensor  of  type 
k  at  node  i,  is  calculated.  Conversely,  the  value  djjk  indicates  the  probability  of  detection  and  can 
be  (but  is  not)  calculated  by  using  djjk  =  1  -  m Also,  given  a  specific  set  of  nodes  i  and  j  in  the 
network,  the  probabilities  and  rrijjk  are  assumed  symmetric  because  the  network  used  in  this 
research  consists  almost  entirely  of  level  farmland,  with  a  uniform  elevation  and  few  obstacles. 
However,  the  probabilities  m///c  and  w7i*  can  differ  in  the  presence  of  obstacles  and  elevation 
differences.  In  order  to  account  for  these  differences,  the  input  parameter  <?//,  ij: 1-673  can  be 
used.  For  networks  with  varying  elevations  and  other  obstacles,  the  values  in  the  e#  matrix  can 
be  set  anywhere  between  0  and  1  allowing  the  probability  of  detection  at  individual  sets  of  i-j 
nodes  to  be  at  its  greatest  value  (e,7=i),  its  lowest  value  (e,7=0),  or  to  be  degraded  ( 0<eij<l ).  The 
level  of  degradation  for  a  particular  i-j  arc  should  be  based  on  the  observed  degradation  created 
by  the  elevation  change  or  obstacle  in  question.  The  formula  for  the  detection  probability  is: 
dijk  =  1  -  mijk  =  (1  -  Beta{5ijt  ak,  pk,  ak,  bfe))e0-;  for  i,j:  1  -  673,  k:  1-4 
But  the  above  formula  is  not  used  in  the  model.  Instead,  the  above  formula  is  used  to 
compute  the  formula  for  miss  probability,  m^,  which  is  then  used  in  the  model. 

mijk  =  1  -  etj  +  Beta(8iJ,ak,pk,ak,blc)eij;  for  i,j:  1  -  673,  k:  1-4 
Furthermore,  since  this  research  assumes  level  ground  and  no  obstacles,  all  values  ejj=l, 
ij: 1-673  for  the  purposes  of  this  research.  With  all  e(/=/, 

rriijk  —  Beta(8ij,  ak,  (3k,  ak,  bkj\  for  i,j:  1  —  673,  k:  1  —  4 
For  combining  probabilities,  it  is  assumed  that  sensor  detections  are  independent,  i.e.  if  a 
sensor  detects  an  intruder  at  a  node  with  probability  a,  and  another  sensor  detects  the  intruder 
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with  a  probability  b,  the  combined  probability  of  detection  =  1-  (1-prob.  of  detect.  a)(l-prob.  of 
detect,  b).  Furthermore,  it  is  assumed  that  some  sensor  types  are  capable  of  both  daytime  and 
nighttime  operation,  while  others  are  capable  of  only  daytime  or  nighttime  operation.  For  each 
node  in  the  network,  probabilities  of  detection  for  daytime  and  nighttime,  df  and  df 
respectively,  are  computed  as  follows: 

df  —  1  —  J~ |  rriijb)  for  i,j:  1  —  673 ,  k\  day 

j,k 

and 

df  =  1  —  m^;  for  i,j:  1  —  673 ,  k\ night 

j,k 

The  daily  probability  of  detection,  /;„  i: 1-673  is  defined  as  the  average  between  the 
daytime  and  nighttime  probabilities  of  detection. 

Pi  —  average{df ,  df );  for  i:  1  —  673 

Average  daily  probability  of  detection,  AvgCov ,  and  minimum  daily  probability  of 
detection,  MinCov,  are  also  defined  as  follows: 

Average  Coverage  —  AvgCov  —  averagefpf) ;  for  i\  1  —  673 
Minimum  Coverage  =  MinCov  =  minimumfpf)  ;  for  i\  1  —  673 

Note  that  the  df  and  df  calculations  are  made  under  the  assumption  of  independent 
probabilities.  Independence  is  a  notional  assumption  used  in  this  and  other  sensor  placement 
literature  (21;  24;  25).  In  actual  practice,  a  certain  amount  of  correlation  may  exist  between 
sensors  (25:4).  However,  if  the  amount  of  correlation  is  determined  to  be  statistically 
insignificant,  the  independent  assumption  can  continue  to  be  used.  Otherwise,  the  formulas  may 
need  to  be  modified  to  account  for  correlation.  Also,  this  research  does  not  address  sensor  fusion, 
i.e.  the  process  by  which  the  data  from  the  various  sensors  is  combined  and  processed. 
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Total  Cost  Calculation. 


The  total  cost  of  building  the  proposed  sensor  network  is  calculated  and  used  as  a 
parameter  forcing  adherence  to  a  budgetary  constraint.  By  design,  sensors  are  placed  at  as  many 
nodes  as  possible  in  order  to  maximize  coverage.  Absent  a  budgetary  constraint,  sensors  would 
be  placed  at  every  node  in  the  network.  There  are  five  cost  parameters.  InfCost  is  the 
infrastructure  cost  for  building  at  a  node.  InfCost  is  the  same  at  every  node  with  sensors  present. 
The  input  parameters  SensCostu,  k:l-4  define  individual  costs  for  purchasing  and  installing  the 
four  sensor  types.  Total  Cost ,  is  the  total  cost  for  building  the  network  and  can  be  calculated 
from  the  five  cost  parameters: 


Total  Cost  =  ^  Si(InfCost  +  ^  sik  *  SensCostk)  ,i\  1  —  673,  k\  1  —  4 


In  the  Total  Cost  calculation,  Si  is  the  binary  selection  variable,  equal  to  1  at  nodes  with 
sensors  and  0  at  nodes  without  sensors,  and  s,k  is  the  binary  input  parameter  describing  which 
sensor  types  are  placed  at  a  node  if  s;=i  at  that  node.  The  total  cost  calculation  is  re-computed  at 
each  iteration. 


Constraints 

The  first  constraint  requires  the  binary  selection  of  nodes  for  sensor  placement,  i.e.  Sj, 
i: 1-67 3  is  binary.  A  number  of  additional  constraints  are  used  to  compel  the  solution  to  exhibit 
desired  attributes. 

Budget. 

A  notional  Budget  is  assumed  to  be  available  for  the  build-up  of  the  network  and  a 
constraint  is  created  so  that  the  total  cost  cannot  exceed  the  total  budget  ( Total  Cost  <  Budget). 
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Node  Availability. 

A  binary  input  parameter  v„  i: 1-673  is  used  to  allow  nodes  to  be  turned  on  or  off  for 
sensor  placement.  The  constraint  .v,-  <  v,  turns  nodes  off  if  v,=0.  The  parameter  v,  is  set  to  1  (on) 
by  default. 

Node  Preference. 

An  input  parameter  w„  i: 1-673  may  be  used  to  require  certain  pre-identified  nodes  to 
have  a  minimum  coverage,  pi.  (25:6)  Note  that  meeting  this  constraint  may  require  more  assets 
than  available  under  the  budget  and  could  result  in  an  infeasible  solution.  For  the  notional 
network  used  to  test  the  model  developed  in  this  research,  there  are  no  preferential  nodes  or 
zones  (sets  of  nodes),  but  the  model  was  designed  to  be  able  to  use  this  constraint. 

Solving  the  Model 

Dhillon  and  Chakrabarty,  describe  two  notional  algorithms  for  sensor  placement  which 
provide  “effective  coverage  and  surveillance  in  distributed  sensor  networks”  (24:1609).  The 
algorithm  called  “MAX_AVG_COV”  is  an  iterative  algorithm  which  places  one  sensor  at  a  time, 
without  backtracking,  until  the  average  miss  probability  drops  below  a  desired  maximum,  i.e.  the 
average  detection  probability  rises  above  a  desired  minimum.  Similarly,  the  algorithm  called 
“MAX_MIN_COV”  is  an  iterative  algorithm  which  places  one  sensor  at  a  time,  without 
backtracking,  until  the  largest  miss  probability  drops  below  a  desired  maximum,  i.e.  the  smallest 
detection  probability  rises  above  a  desired  minimum. 

The  model  developed  in  this  research  is  solved  using  algorithms  similar  to  the  ones 
described  by  Dhillon  and  Chakrabarty.  Additionally,  the  Premium  Solver  Platform  for  Excel, 
together  with  the  OptQuest  field-installable  engine,  is  used  to  solve  the  model. 
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VBA  (Visual  Basic  for  Applications) 

Two  algorithms,  VBA-AvgCov  and  VBA-MinCov,  were  developed  using  the  Visual  Basic 
for  Applications  (VBA)  programming  language  in  Excel.  VBA-AvgCov  is  an  iterative  greedy 
algorithm  which  places  one  sensor  package  at  a  time  on  the  network  (selects  one  node  at  a  time) 
until  all  funds  are  exhausted,  i.e.  when  the  Total  Cost  exceeds  the  Budget ,  the  algorithm  stops 
and  returns  the  previous  (last  feasible)  solution.  At  each  iteration,  the  algorithm  places  a  sensor 
package  at  the  node  that  will  affect  the  greatest  incremental  increase  in  the  average  coverage, 
AvgCov.  There  is  no  backtracking  in  this  algorithm.  VBA-MinCov  works  in  a  similar  fashion  but 
always  chooses  the  node  that  will  produce  the  greatest  incremental  increase  in  the  minimum,  as 
opposed  to  the  average,  coverage  MinCov.  However,  since  the  first  few  iterations  are  likely  to 
produce  MinCov=0,  because  there  are  not  yet  enough  sensors  to  cover  all  nodes,  the  algorithm 
chooses  a  sensor  location  which  maximizes  AvgCov  until  it  reaches  a  point  where  MinCov>0. 
Then,  for  the  first  iteration  where  MinCov>0 ,  and  each  subsequent  iteration,  the  algorithm 
switches  to  choosing  a  sensor  location  at  the  node  which  maximizes  MinCov.  This  algorithm  also 
works  without  backtracking.  The  VBA  code  for  both  algorithms  can  be  found  in  Appendix  B. 

OptQuest  Solver. 

The  two  algorithms  described  above  are  used  as  a  baseline  for  the  model  However,  in  an 
effort  to  ensure  good  results,  an  attempt  is  made  to  improve  upon  the  baseline  solutions  using 
commercial  solver  software.  Unfortunately,  due  to  its  size  and  complexity,  the  model  cannot  be 
solved  using  the  built-in  Excel  Solver.  In  fact,  the  model  exceeds  the  number  of  variables  and 
constraints  that  the  Excel  Solver  can  handle.  It  is  also  non-linear  and  non-smooth  (due  mostly  to 
the  use  of  ‘min’  and  ‘if  functions)  and  the  Excel  Solver  requires  linearity. 
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However,  OptQuest  solver  engine  from  Frontline  Systems,  the  makers  of  the  Excel 

Solver,  “employs  metaheuristics  such  as  tabu  search  and  scatter  search  to  solve  nonsmooth 

optimization  problems  of  up  to  5,000  variables  and  1,000  constraints.”  (27) 

In  the  tabu  search  category  of  meta-heuristics,  the  essential  idea  is  to  'forbid'  search 
moves  to  points  already  visited  in  the  (usually  discrete)  search  space,  at  least  for  the 
upcoming  few  steps.  That  is,  one  can  temporarily  accept  new  inferior  solutions,  in  order 
to  avoid  paths  already  investigated.  This  approach  can  lead  to  exploring  new  regions  of  D 
[the  search  space],  with  the  goal  of  finding  a  solution  by  'globalized'  search.  Tabu  search 
has  traditionally  been  applied  to  combinatorial  optimization  (e.g.,  scheduling,  routing, 
traveling  salesman)  problems.  (33) 

Scatter  search  operates  on  a  set  of  solutions,  the  reference  set,  by  combining  these 
solutions  to  create  new  ones.  The  main  mechanism  for  combining  solutions  is  such  that  a 
new  solution  is  created  from  the  linear  combination  of  two  other  solutions.  The  reference 
set  may  evolve  [over  time].  (34) 

In  order  to  use  the  OptQuest  solver  engine,  two  software  packages  need  to  be  installed: 
the  first  is  Frontline  Systems’  Premium  Solver  Platform  (PSP),  and  the  second  is  the  OptQuest 
solver  engine  itself.  The  combined  package  is  simply  referred  to  as  the  “OptQuest  solver.”  While 
the  OptQuest  solver  can  only  guarantee  optimality  by  complete  enumeration,  Frontline  Systems 
claims  that  the  OptQuest  solver  finds  “remarkably  good  solutions  with  unprecedented  speed” 
(27).  In  Chapter  4,  this  claim  is  evaluated  by  comparing  the  solver  results  against  the  results  of 
the  VBA  algorithms. 
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Review 


The  models  discussed  above  can  be  summarized  as  follows: 

Max  AvgCov  ( Maximize  Average  Coverage ) 

s.t.  subject  to, 

Total  Cost  <  Budget  (budgetary  constraint) 

st  <  vt  (node  on/ojf  constraint) 

wt  <  Pi  (minimum  coverage  constraint) 

And 

Max  MinCov  (Maximize  Minimum  Coverage) 

s.t.  subject  to, 

Total  Cost  <  Budget  (budgetary  constraint) 

st  <  Vi  (node  on/ojf  constraint) 

wt  <  Pi  (minimum  coverage  constraint) 

There  is  one  model  with  two  objectives.  The  VBA  algorithms  developed  above,  and  the 
OptQuest  Solver,  are  used  to  solve  the  optimization  model  for  each  objective.  AvgCov  represents 
the  average  daily  coverage  over  the  entire  network  and  MinCov  represents  the  node  in  the 
network  which  has  the  least  sensor  coverage.  The  Total  Cost  is  calculated  by  adding  up 
infrastructure  and  sensor  costs  for  all  selected  nodes  in  the  network.  The  Budget  is  an  input 
parameter.  For  nodes  i  =  1-673,  Si  is  a  binary  variable  indicating  selected  nodes,  v,  is  a  binary 
input  parameter  which,  when  set  equal  to  0,  prevents  node  selection,  pi  is  the  calculated  value  for 
daily  coverage  (as  defined  previously),  and  w;  is  an  input  parameter  used  to  force  minimum  daily 
coverage. 
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IV.  Results  and  Conclusions 


In  Chapter  3,  a  model  is  developed  for  placement  of  distributed  sensors  on  a  network 
with  the  goals  of  maximizing  minimum  coverage  and  average  coverage  on  the  network.  In 
Chapter  4,  the  model  is  tested  using  the  methods  proposed  in  Chapter  3.  First,  the  VBA 
algorithms  are  tested.  Then,  the  OptQuest  solver  is  tested  and  the  results  are  compared  to  the 
previous  results.  The  computer  used  for  the  tests  is  an  Asus  A8jp  laptop  with  an  Intel  Core  2  Duo 
7200  processor  and  2GB  of  RAM  running  Microsoft  Excel  2007  and  OptQuest  solver  7.0.  Run 
times  are  preserved  for  each  of  the  runs  but  the  computer  must  be  used  for  other  work  at  the 
same  time  as  the  Excel  runs.  Therefore,  the  quoted  times  may  not  represent  the  full  capability  of 
the  computer  being  used. 

Inputs 

A  number  of  inputs  are  required  to  run  the  model.  The  673  node  network  developed  in 
Chapter  3  is  the  primary  data  input  to  the  model.  Figure  10  shows  a  representation  of  this 
network  in  Google  Earth  with  circles  representing  intersection  nodes  and  diamonds  representing 
centers  of  land  parcels. 
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Figure  10.  Complete  Network 


37 


However,  the  full  network,  shown  in  Figure  10,  contains  a  total  of  47  nodes  which  are  not 
included  in  the  optimization  for  one  of  two  reasons.  A  set  of  4  nodes  (orange  nodes  in  Figure  11) 
are  excluded  because  they  are  located  away  from  the  rest  of  the  network,  and  a  set  of  43  nodes 
(red  nodes  in  Figure  1 1 )  are  excluded  because  they  are  located  in  the  city.  It  is  assumed  that 
sensors  will  not  be  effective  inside  a  city.  For  the  i  values  corresponding  to  the  eliminated  nodes, 
the  parameter  v,  is  set  to  0  (vi  =  0).  The  remaining  626  nodes  have  v,-  =  1. 


Figure  11.  Network  with  eliminated  nodes  highlighted 


Parameter  Inputs. 

The  parameters  chosen  for  these  optimization  runs  (see  Table  1  and  Table  2)  cannot  be 
validated  against  true  operational  settings.  Instead,  the  notional  parameters  selected  appear 
“reasonable”  for  the  purpose  of  these  tests.  For  example,  the  4  kilometer  range  assumption  for 
the  Ground  Search  Radar  (sensor  3)  is  in  line  with  the  available  data  (3).  The  Budget  parameter 
is  varied  in  order  to  test  the  model  under  differing  conditions. 
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Table  1.  Sensor  Input  Parameters 


Sensor 

Sensor 

Day/Night 

Shape  Parameters 

Bounds 

Cost 

Name 

k 

day  night 

ak 

Pk 

ak  bk 

SensCostk 

Regular  Camera 

HI 

yes  no 

0  3 

$5000 

Infrared  Camera 

□ 

no  yes 

mi 

0  1.2 

$4000 

Ground  Search  Radar 

a 

yes  yes 

1.2 

0.8 

0  4 

$6000 

Seismic  Sensor 

D 

yes  yes 

0.6 

1.3 

0  0.8 

$2000 

Table  2.  Additional  Inputs 


Total  Construction  Budget 

Budget 

$100,000;  $300,000 

Infrastructure  Cost  Per  Node 

InfCost 

$7000 

It  is  assumed  that  seismic  sensors  ( k=4 )  cannot  function  effectively  when  placed  at 
intersection  nodes.  For  this  reason,  a  seismic  sensor  will  not  be  part  of  the  sensor  package  if  an 
intersection  node  is  selected.  Table  3  shows  that  the  sensor  types  1,  2,  and  3  (su,  sq,  and  su)  will 
be  placed  at  intersection  nodes  (i:  i= 1-399)  while  all  4  sensor  types  will  be  placed  at  center  of 
land  parcel  nodes  (i:  i=400-673). 


Table  3.  Available  Sensors 


Node  Type 

i 

Sil 

Si2 

Si3 

Si4 

Intersection  Node 

1-399 

1 

1 

1 

0 

Center  of  Land  Parcel 

400-673 

1 

1 

1 

1 

Lastly,  due  to  the  problem  assumptions,  as  explained  in  Chapter  3,  the  parameters  ey  and 
Wi  where  ij: 1-673,  were  set  to  one  and  zero  respectively. 

Detection  Curves. 

The  shape  and  bound  parameters  from  Table  1  produce  the  following  detection  curves  for 
the  four  sensor  types: 
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Sensor  1  Detection  Curve 


Figure  12.  Sensor  1  Detection  Probability  Curve 


Figure  13.  Sensor  2  Detection  Probability  Curve 
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Sensor  3  Detection  Curve 
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Figure  14.  Sensor  3  Detection  Probability  Curve 


Figure  15.  Sensor  4  Detection  Probability  Curve 

Figures  12  through  15  show  the  detection  and  miss  curves  for  each  of  the  four  sensor 
types,  over  their  individual  ranges.  The  curves  are  used  to  calculate  point  to  point  miss 
probabilities  for  each  pair  of  nodes  in  the  network  and  for  each  sensor  type. 
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Results 


Every  test  is  completed  two  times;  once  assuming  a  budget  of  $100,000  {Run  1),  and  a 
second  time  assuming  a  budget  of  $300,000  ( Run  2).  The  results  of  each  of  these  two  runs  is 
shown  below: 

VBA-AvgCov  Algorithm  -  Run  1. 

Table  4.  VBA-AvgCov  Algorithm  Run  1  Summary 


Total  Construction  Budget 

Budget 

$100,000 

Total  Construction  Cost 

Total  Cost 

$96,000 

Number  of  Sensors  Selected 

Ysi 

4 

Objective  (Maximization) 

Average(pi) 

0.7200 

Run  Time 

Hours 

3  (approx.) 

Table  5.  VBA-AvgCov  Algorithm  Run  1  Probability  Summary 


Sensor 

Day/Night 

Variable 

Average 

(Max)/Min 

Std.  Dev. 

Probability  of  Missing 

1 

Day 

mu 

0.4268 

(1.0000) 

0.3280 

2 

Night 

mi2 

0.9547 

(1.0000) 

0.1378 

3 

Both 

mi3 

0.3420 

(1.0000) 

0.2410 

4 

Both 

mi4 

0.9826 

(1.0000) 

0.0935 

Probability  of  Detection 

1,2,4 

Day 

did 

0.7795 

0.0000 

0.2668 

2,3,4 

Night 

din 

0.6605 

0.0000 

0.2439 

1,2, 3, 4 

Average 

Pi 

0.7200 

0.0000 

0.2525 
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Figure  16.  VBA-AvgCov  Algorithm  Run  1  Visual  Summary 


Figure  16  shows  the  4  nodes  selected  in  Run  1  highlighted  in  green.  The  selected  sensors 
are  placed  in  a  relatively  flat  line  across  the  center  of  the  network.  This  selection  produces  an 
average  coverage  across  the  network  of  72%.  However,  there  are  12  nodes,  highlighted  in 
yellow,  which  are  not  covered  by  any  of  the  sensors.  Due  to  their  location,  away  from  the  border 
and  at  the  very  back  of  the  coverage  area,  these  uncovered  nodes,  while  not  ideal,  are  not  a  great 
cause  for  concern. 

VBA-AvgCov  Algorithm  -  Run  2. 

Table  6.  VBA-AvgCov  Algorithm  Run  2  Summary 


Total  Construction  Budget 

Budget 

$300,000 

Total  Construction  Cost 

Total  Cost 

$286,000 

Number  of  Sensors  Selected 

5>/ 

12 

Objective  (Maximization) 

Average(pi) 

0.9700 

Run  Time 

Hours 

8  (approx.) 
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Table  7.  VBA-AvgCov  Algorithm  Run  2  Probability  Summary 


Sensor 

Day/Night 

Variable 

Average 

(Max)/Min 

Std.  Dev. 

Probability  of  Missing 

1 

Day 

mu 

0.0870 

(1.0000) 

0.1286 

2 

Night 

mi2 

0.8721 

(1.0000) 

0.2117 

3 

Both 

m3 

0.0505 

(0.3980) 

0.0598 

4 

Both 

mi4 

0.9519 

(1.0000) 

0.1537 

Probability  of  Detection 

1,2,4 

Day 

did 

0.9888 

0.7183 

0.0316 

2,3,4 

Night 

din 

0.9513 

0.6020 

0.0603 

1,2, 3, 4 

Average 

Pi 

0.9700 

0.6643 

0.0445 

Figure  17.  VBA-AvgCov  Algorithm  Run  2  Visual  Summary 


Since  the  VBA-AvgCov  algorithm  does  not  backtrack,  both  runs  are  identical  until  the 
fourth  sensor.  However,  while  Run  1  ended  after  placing  the  fourth  sensor  due  to  the  smaller 
budgetary  constraint,  in  Run  2  the  algorithm  was  allowed  to  continue  placing  sensors  until  the 
larger  budget  was  exhausted.  The  Run  2  solution  improves  the  average  coverage  by  25%  over  the 
Run  1  solution.  Also,  there  are  no  longer  any  nodes  with  zero  coverage.  The  node  shown  in 
yellow  in  Figure  17  is  the  node  with  the  lowest  average  coverage  at  0.664. 
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VBA-MinCov  Algorithm  -  Run  1. 

The  VBA-MinCov  algorithm  works  the  same  way  as  the  VBA-AvgCov  algorithm  until 
there  is  a  minimum  probability  of  detection  greater  than  zero;  at  which  point  the  algorithm  starts 
to  choose  sensor  locations  in  order  to  maximize  the  minimum  probability  of  detection  at  each 
iteration.  At  the  $100,000  budget  level  this  threshold  is  not  met,  so  Run  1  for  the  VBA-MinCov 
algorithm  is  identical  to  Run  1  for  the  VBA-AvgCov  algorithm.  Run  2,  however,  does  produce 
results  that  are  significantly  different  from  those  produced  by  Run  2  of  the  VBA-AvgCov 
algorithm. 

VBA-MinCov  Algorithm  -  Run  2. 

Table  8.  VBA-MinCov  Algorithm  Run  2  Summary 


Total  Construction  Budget 

Budget 

$300,000 

Total  Construction  Cost 

Total  Cost 

$294,000 

Number  of  Sensors  Selected 

I  Si 

13 

Objective  (Maximization) 

Minimum(pi) 

.5156 

Run  Time 

Hours 

8  (approx.) 

Table  9.  VBA-MinCov  Algorithm  Run  2  Probability  Summary 


Sensor 

Day/Night 

Variable 

Average 

(Max)/Min 

Std.  Dev. 

Probability  of  Missing 

1 

Day 

mu 

0.1303 

(1.0000) 

0.1905 

2 

Night 

ml2 

0.8827 

(1.0000) 

0.2159 

3 

Both 

mi3 

0.0831 

(0.5684) 

0.1098 

4 

Both 

mi4 

0.9526 

(1.0000) 

0.1595 

Probability  of  Detection 

1,2,4 

Day 

did 

0.9698 

0.5179 

0.0722 

2,3,4 

Night 

din 

0.9183 

0.4316 

0.1104 

1,2, 3, 4 

Average 

Pi 

0.9441 

0.5156 

0.0893 
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Figure  18.  VBA-MinCov  Algorithm  Run  2  Visual  Summary 


Table  9  shows  the  numerical,  and  Figure  18  shows  the  visual,  result  of  the  VBA-MinCov 
algorithm.  Both  numerically  and  visually,  the  results  of  the  VBA-MinCov  algorithm  are  inferior 
to  those  provided  by  the  VBA-AvgCov  algorithm.  Using  the  VBA-MinCov  algorithm,  the  average 
coverage  is  approximately  2.5%  lower  while  minimum  coverage  is  approximately  15%  lower 
than  the  results  gained  from  the  VBA-AvgCov  algorithm.  The  yellow  node  in  Figure  18  is  the 
node  with  the  lowest  average  coverage  (pi)  at  0.516.  These  results  imply  that  the  VBA-AvgCov 
algorithm  is  superior  to  the  VBA-MinCov  algorithm  even  when  the  objective  is  to  maximize 
minimum  coverage. 

OptQuest  Solver 

The  OptQuest  solver  (OQS)  for  Excel  is  used  in  an  attempt  to  improve  upon  the  results  of 
the  two  VBA  algorithms  seen  above.  The  OptQuest  solver  is  given  the  same  objectives  and 
constraints  as  the  VBA  algorithms  and  two  runs  are  completed  for  each  objective  just  as  before. 
All  runs  are  completed  using  the  OptQuest  solver  with  30,000  iterations  unless  mentioned 
otherwise.  Appendix  C  shows  all  of  the  OptQuest  solver  inputs  and  options  used  for  the 
following  runs. 
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OQS-AvgCov  -  Run  1. 


Table  10.  OQS-AvgCov  Run  1  Summary 


Total  Construction  Budget 

Budget 

$100,000 

Total  Construction  Cost 

Total  Cost 

$90,000 

Number  of  Sensors  Selected 

Yst 

4 

Objective  (Maximization) 

Average(pi) 

0.7403 

Run  Time 

Hours 

11:04:31 

Table  11.  OQS-AvgCov  Run  1  Probability  Summary 


Sensor 

Day/Night 

Variable 

Average 

(Max)/Min 

Std.  Dev. 

Probability  of  Missing 

1 

Day 

mu 

0.4064 

(1.0000) 

0.3049 

2 

Night 

mi2 

0.9511 

(1.0000) 

0.1421 

3 

Both 

mi3 

0.3303 

(1.0000) 

0.2053 

4 

Both 

mi4 

0.9813 

(1.0000) 

0.0951 

Probability  of  Detection 

1,2,4 

Day 

did 

0.8077 

0.0000 

0.2327 

2,3,4 

Night 

din 

0.6729 

0.0000 

0.2092 

1,2, 3, 4 

Average 

Pi 

0.7403 

0.0000 

0.2175 

Figure  19.  OQS-AvgCov  Run  1  Visual  Summary 


After  30,000  iterations,  the  OptQuest  solver  found  a  solution  with  0.740  average 
coverage;  approximately  a  2%  improvement  over  the  VBA-AvgCov  Run  1  result.  Also,  while  the 


47 


OptQuest  solver  did  not  find  a  solution  with  minimum  coverage  greater  than  zero,  it  did  reduce 
the  number  of  uncovered  nodes  (yellow  nodes  in  Figure  19)  to  2,  down  from  12  under  the  VBA- 
AvgCov  Rim  1  result.  The  standard  deviation  is  also  lower  using  the  OptQuest  solver. 

OQS-AvgCov  -  Run  2. 


Table  12.  OQS-AvgCov  Run  2  Summary 


Total  Construction  Budget 

Budget 

$300,000 

Total  Construction  Cost 

Total  Cost 

$298,000 

Number  of  Sensors  Selected 

Yst 

13 

Objective  (Maximization) 

Average(pi) 

0.9795 

Run  Time 

Hours 

9:58:32 

Table  13.  OQS-AvgCov  Run  2  Probability  Summary 


Sensor 

Day/Night 

Variable 

Average 

(Max)/Min 

Std.  Dev. 

Probability  of  Missing 

1 

Day 

mu 

0.0574 

(0.4269) 

0.0618 

2 

Night 

m  12 

0.8436 

(1.0000) 

0.2269 

3 

Both 

mu 

0.0387 

(0.3264) 

0.0390 

4 

Both 

mi4 

0.9376 

(1.0000) 

0.1675 

Probability  of  Detection 

1,2,4 

Day 

did 

0.9957 

0.9012 

0.0095 

2,3,4 

Night 

din 

0.9638 

0.6736 

0.0395 

1,2, 3, 4 

Average 

Pi 

0.9797 

0.7874 

0.0241 
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Figure  20.  OQS-AvgCov  Run  2  Visual  Summary 


The  OptQuest  solver  produced  a  solution  with  a  1.5%  increase  in  average  coverage,  and 
12%  increase  in  minimum  coverage,  over  the  VBA-AvgCov  algorithm  Run  2  solution.  The  yellow 
node  in  Figure  20  is  the  node  with  the  lowest  average  coverage  (pi)  at  0.787. 

OQS-MinCov  -  Run  1. 

After  30,000  iterations,  the  OptQuest  solver  was  not  able  to  find  a  solution  with 
minimum  coverage  greater  than  zero  (given  the  $100,000  budgetary  constraint  for  Run  1 ).  The 
OptQuest  solver  can  only  search  for  one  objective  at  a  time  and,  failing  to  find  a  solution  with 
MinCov  >  0,  it  simply  returned  an  all  zero  solution  set,  i.e.  no  nodes  were  selected  for  sensor 
placement.  While  this  failure  cannot  guarantee  that  a  solution  with  MinCov  >  0  does  not  exist,  it 
is  part  of  a  larger  pattern  which  appears  to  point  in  that  direction.  However,  the  OptQuest  solver 
can  only  prove  MinCov  =  0  by  complete  enumeration.  Therefore,  additional  runs  may  be  able  to 
find  a  solution  with  MinCov  >  0. 
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OQS-MinCov  -  Run  2. 


Table  14.  OQS-MinCov  Run  2  Summary 


Total  Construction  Budget 

Budget 

$300,000 

Total  Construction  Cost 

Total  Cost 

$294,000 

Number  of  Sensors  Selected 

13 

Objective  (Maximization) 

Minimum(pi) 

0.8244 

Run  Time 

Hours 

10:53:11 

Table  15.  OQS-MinCov  Run  2  Probability  Summary 


Sensor 

Day/Night 

Variable 

Average 

(Max)/Min 

Std.  Dev. 

Probability  of  Missing 

1 

Day 

mu 

0.0679 

(0.5777) 

0.0750 

2 

Night 

mi2 

0.8417 

(1.0000) 

0.2290 

3 

Both 

mi3 

0.0437 

(0.2404) 

0.0397 

4 

Both 

mi4 

0.9368 

(1.0000) 

0.1668 

Probability  of  Detection 

1,2,4 

Day 

did 

0.9944 

0.8892 

0.0116 

2,3,4 

Night 

din 

0.9590 

0.7596 

0.0406 

1,2, 3, 4 

Average 

Pi 

0.9767 

0.8244 

0.0255 

Figure  21.  OQS-MinCov  Run  2  Visual  Summary 


The  OQS-MinCov  Run  2  solution  has  a  3.3%  higher  average  coverage  and  a  30.9% 
higher  minimum  coverage  than  the  VBA-MinCov  algorithm  Run  2  solution.  Also,  the  OQS- 
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MinCov  Run  2  solution  has  an  average  coverage  almost  identical  to  the  OQS-AvgCov  Run  2 
solution  (0.3%  smaller)  and  a  3.7%  higher  minimum  coverage  than  the  OQS-AvgCov  Run  2 
solution.  The  yellow  node  in  Figure  21  is  the  node  with  the  lowest  average  coverage  (pi)  at 
0.824. 

Summary 

The  OptQuest  solver  produced  consistently  better  results  than  the  VBA  algorithms.  Table 
16  shows  that,  under  the  Max  AvgCov  objective,  the  OptQuest  solver  produced  the  highest 
average  coverage  and,  under  the  Max  MinCov  objective,  the  OptQuest  solver  produced  the 
highest  minimum  coverage.  The  VBA  algorithms  produced  worse  results  in  both  cathegories. 


Table  16.  Result  Summary 


Objective 

Type 

AverageQ?,) 

Minimum^,) 

Std.  Dev.(pf) 

Run  1 

Max-Avg 

VBA 

0.7200 

0.0000 

0.2525 

Solver 

0.7403 

0.0000 

0.2175 

Max-Min 

VBA 

0.7200 

0.0000 

0.2525 

Solver 

0.0000 

0.0000 

0.0000 

Run  2 

Max-Avg 

VBA 

0.9700 

0.6643 

0.0445 

Solver 

0.9797 

0.7874 

0.0241 

Max-Min 

VBA 

0.9441 

0.5156 

0.0893 

Solver 

0.9767 

0.8244 

0.0255 

Model  Behavior  Analysis 
Feasibility. 

Under  the  $100,000  budgetary  assumption,  the  OptQuest  solver  cannot  find  a  feasible 
solution  for  the  maximization  of  the  minimum  coverage  ( OQS-MinCov  Run  1 ).  The  budget 
allows  only  4  sensor  packages  to  be  purchased  and,  with  4  sensor  locations,  the  objective,  Max 
MinCov,  appears  to  equal  zero,  i.e.  there  is  always  at  least  one  node  which  is  not  covered  by  the 
solution.  Since  it  cannot  find  a  feasible  solution,  the  OptQuest  solver  returns  an  empty  solution 
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set,  i.e.  no  nodes  are  selected  for  sensor  placement.  Furthermore,  although  the  VBA-MinCov 
algorithm  does  produce  an  answer  under  the  $100,000  budgetary  assumption  (Run  1),  the  answer 
is  based  solely  on  the  secondary  objective  (Max  AvgCov)  and  is  identical  to  the  answer  produced 
by  Run  1  of  the  VBA-AvgCov  algorithm. 

Node  Separation. 

Under  the  $100,000  budgetary  assumption,  sensors  are  approximately  equidistant  from 
each  other,  while  under  the  $300,000  budgetary  assumption,  sensors  are  placed  both  close 
together  and  far  apart.  This  leads  to  a  visual  solution,  Figure  21,  which  appears  to  have  large 
gaps  in  it.  However,  due  to  overlapping  coverage  created  by  additional  sensors,  the  apparent 
gaps  are  actually  covered  very  well;  a  fact  which  is  attested  to  by  a  Standard  Deviation  value  of 
approximately  2.5%.  Figure  22  shows  the  coverage  for  sensor  4  (Ground-Search  Radar)  for  the 
OQS-MinCov  Run  2  solution.  The  apparent  gaps  in  Figure  21  are  very  well  explained  by  the 
coverage  map  shown  in  Figure  22.  The  figure  also  explains  the  low  standard  deviation. 


Figure  22.  OQS-MinCov  Run  2  Sensor  4  Coverage  Map 
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Figure  22  shows  that  the  number  of  sensors  covering  a  particular  node  is  directly  related  to  that 


node’s  distance  from  the  closest  sensor.  Since  a  sensor’s  probability  of  detection  drops  as  the 
distance  from  the  sensor  increases,  the  solution  covers  the  nodes  that  are  farther  from  a  sensor 
with  additional  sensors. 

Solution  Speed. 

The  four  OptQuest  solver-based  solutions  shown  above  took  30,000  iterations,  and  as 
long  as  1 1  hours,  to  obtain.  In  some  instances,  it  may  be  necessary  to  obtain  solutions  faster  than 
they  can  be  obtained  with  30,000  iterations.  Four  additional  runs  of  the  OptQuest  solver  with  the 
Max  AvgCov  objective  in  order  to  test  the  solver’s  ability  to  find  good  solutions  quickly.  The 
additional  runs  were  set  to  500,  1000,  5000,  and  10,000  iterations  respectively.  Table  17  shows 
the  results  of  these  runs  as  well  as  the  original  result  for  30,000  iterations.  A  budgetary  constraint 
of  $300,000  was  used  for  all  runs. 

Table  17.  OQS-AvgCov  Iteration  Comparison 


Iterations 

AvgCov 

MinCov 

Std.  Dev. 

Run  Time 

500 

0.9648 

0.5870 

0.0620 

00:34:21 

1000 

0.9735 

0.6737 

0.0426 

00:36:40 

5000 

0.9709 

0.6314 

0.0435 

01:49:17 

10000 

0.9787 

0.8026 

0.0270 

04:53:02 

30000 

0.9797 

0.7874 

0.0241 

09:58:32 

Table  17  shows  that  a  good  result  can  be  achieved  within  as  few  as  1000  iterations.  Also 
note  that  the  result  after  1000  iterations  is  better  than  the  result  after  5000  iterations.  Since  the 
OptQuest  solver  is  allowed  to  pick  different  random  number  sequences  with  each  run  (see 
Appendix  C,  Figure  28)  the  results  can  and  do  differ  with  each  run.  However,  the  general  trend 
shows  that  the  results  improve  with  the  number  of  iterations.  If  time  is  not  of  the  essence,  it  is 
suggested  that  at  least  30,000  iterations  be  completed. 
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Conclusion 


The  research  developed  a  “proof-of-concept”  model  for  distributed  sensor  placement 
optimization  for  border  security  using  Microsoft  Excel  and  Frontline  Systems’  OptQuest  solver. 
The  road  network  of  Calexico,  California,  which  was  used  for  this  model,  was  manually  created 
using  Google  Inc.’s  Google  Earth  application  and  transferred  to  Excel  using  the  free  GEpath 
application.  The  model  optimizes  the  placement  of  electronic  sensors,  in  order  to  maximize  the 
average  per-arc  probability  of  detection  over  a  network,  given  a  budgetary  constraint. 

The  model  was  tested  using  two  separate  objectives;  the  first  is  the  Max  AvgCov 
objective,  which  maximizes  the  average  probability  of  detection  over  all  nodes  of  the  network, 
and  the  second  is  the  Max  MinCov  objective,  which  maximizes  the  minimum  probability  of 
detection  at  any  node  in  the  network.  Each  of  the  two  objectives  was  tested  using  Visual  Basic 
for  Applications  (VBA)  algorithms  and  Frontline  Systems’  Premium  Solver  Platform  with  the 
OpQuest  solver  engine  (OQS). 

Of  the  two  VBA  algorithms,  VBA-AvgCov  and  VBA-MinCov ,  the  VBA-AvgCov  algorithm 
produced  much  better  results.  However  the  OptQuest  solver  was  able  to  produce  quicker  and 
better  results  than  both  of  the  VBA  algorithms. 

Using  the  OptQuest  solver,  the  Max  MinCov  objective  produced  good  results  under  the 
larger  budgetary  constraint,  but  could  not  find  a  feasible  solution  under  the  smaller  budgetary 
constraint.  The  Max  AvgCov  objective  produced  good  results  regardless  of  the  budgetary 
constraint.  Both  objectives  resulted  in  good  sensor  placement  solutions  for  coverage 
maximization  and  are  recommended  for  further  research. 
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Military  Application 

Due  to  the  ease  of  use  of  the  model  (it  runs  on  Microsoft  Excel)  and  its  portability  (it  can 
run  on  any  laptop  with  the  needed  software)  the  model  could  be  used  by  military  forces  in  the 
field  to  determine  positioning  of  sensors  for  border  and  perimeter  security.  While  Google  Earth 
was  used  to  create  the  network  for  this  model,  troops  in  the  field  may  not  have  access  to  this 
software  due  to  the  lack  of  internet  connectivity.  Nonetheless,  the  network  can  be  just  as  easily 
created  with  any  software  capable  of  displaying  latitude  and  longitude  coordinates. 
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V.  Recommendations 


Sensor  Locations  and  Ranges 

The  model  assumes  multiple  sensor  types  are  placed  together  as  packages  in  one  location. 
This  assumption  creates  fewer  physical  sensor  locations  making  it  easier  and  cheaper  to  secure 
the  network  than  if  each  sensor  is  allowed  its  own  location.  However,  each  sensor  has  a  different 
range  and,  under  the  stated  assumption,  the  sensors  with  the  longest  range  overwhelmingly 
influences  the  placement  of  sensors,  i.e.  the  longer  the  range  of  a  sensor,  the  more  nodes  it 
detects.  So,  in  effect  optimization  is  prioritized  based  on  the  range  of  each  sensor,  with  the 
longest  range  having  the  greatest  optimization  priority. 

In  order  to  improve  upon  the  present  solution,  each  sensor  type  will  need  to  be  selected, 
at  each  node,  independent  of  the  other  sensor  types.  This  will  require  either  a  quadrupling  of  the 
number  of  variables,  constraints,  and  calculations  in  the  model,  or  a  separate  budgetary 
constraint  for  each  sensor  type.  If  each  sensor  type  is  given  its  own  budgetary  constraint,  a 
simplified  version  of  this  model  can  be  solved  separately  for  each  sensor  type.  This  will  likely 
result  in  many  more  nodes  with  at  least  one  type  of  sensor  and  the  cost  of  building  the  network 
infrastructure  may  increase  significantly. 

Sensor  Fusion 

The  model  assumes  independence  and  combines  all  probabilities  using  this  assumption. 
However,  in  an  operational  environment,  the  assumption  of  independence  may  not  be  justifiable. 
There  may  be  correlation  between  sensors  of  different  types  as  well  as  between  sensors  of  the 
same  type  at  different  locations.  There  may  also  be  a  number  of  additional  issues  regarding  the 
fusion  of  data  between  individual  sensors,  which  will  need  to  be  considered  and  integrated  into 
the  model  architecture. 
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Summary 

This  research  has  developed  a  proof-of-concept  sensor  placement  model  for  border 
interdiction.  The  model  allows  each  node  to  easily  be  turned  on  or  off  for  sensor  placement 
based  on  economic,  political,  and  operational  considerations.  It  also  has  the  capability  to 
degrade,  or  turn  off,  specific  node  to  node  interactions.  This  capability  allows  more  accurate 
models  of  locations  with  many  obstacles  and/or  elevation  changes  to  be  created.  However,  the 
model  has  not  been  tested  with  operational  inputs  and  needs  to  be  modified  before  it  can  be 
employed  in  an  operational  environment. 
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Appendix  A 


Beta  Probability  Distribution 

The  “beta  [probability]  distribution  is  a  two-parameter  family  of  continuous  probability 
distributions  defined  on  the  interval  [0,  1]”  (8).  By  changing  the  parameters,  a  and  p,  the  Beta 
distribution  can  exhibit  an  infinite  number  of  density  function  shapes.  (8;  35;  36:178,179) 


/  O;  a,  /?) 


Beta  Probability  Density  Function  (pdf): 

xa-\l-x)P-i  _r(a  +  P)  a_1 


■  x 


JgUa  1  (1  —  u)P  1du  r(a)r(P) 

T  is  the  gamma  function  and  B  is  the  beta  function 
Beta  Density  Function  Shapes: 

•  a  <  1,  P  <  1  is  U-shaped 

•  a<l,P>lora=l,P>lis  strictly  decreasing 

o  a  =  1 ,  P  >  1  is  strictly  convex 
o  a=l,P  =  2isa  straight  line 
o  a=l,  l<P<2is  strictly  concave 

•  a  =  1 ,  P  =  2  is  the  uniform  distribution 

•  a=l,P<lora>l,P<lis  strictly  increasing 

o  a  >  2,  P  =  1  is  strictly  convex 
o  a  =  2,  P  =  1  is  a  straight  line 
o  l<a<2,  P  =  1  is  strictly  concave 

•  a  >  1 ,  P  >  1  is  unimodal 

•  If  a  =  p  then  the  density  function  is  symmetric  about  Vi 


B(a,/?) 


x 


a-1(l  —  x)^-1 
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Examples 


Figure  23.  Probability  density  functions  (8) 


Figure  24.  Cumulative  distribution  function  (8) 
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Appendix  B 


The  Visual  Basic  for  Applications  (VBA)  code  for  the  Max_Ave_Det  and  VBA-MinCov 
algorithms  is  presented  below: 

VBA-AvgCov  Algorithm 

Sub  MyAveSolver() 

Dim  MyCurNode  As  Range 

Set  MyCurNode  =  Worksheets("IO").Range("T31") 

Dim  MyFirstCell  As  Range 

Set  MyFirstCell  =  Worksheets("IO").Range("E2") 

Dim  MyLastCell  As  Range 

Set  MyLastCell  =  Worksheets("IO").Range("E674") 

Dim  MyTarCell  As  Range 

Set  MyTarCell  =  Worksheets("IO").Range("T19") 

Dim  MyBudget  As  Range 

Set  MyBudget  =  Worksheets("IO").Range("N13") 

Dim  MyCost  As  Range 

Set  MyCost  =  Worksheets("IO").Range("N23") 

Dim  MyPrev  As  Range 

Set  MyPrev  =  Worksheets("IO")-Range("T29") 

Dim  MyNext  As  Range 

Set  MyNext  =  Worksheets("IO")-Range("T30") 

Dim  MyTotal  As  Integer 
MyTotal  =  0 

Application. Calculation  =  xlManual 
Do 

MyFirstCell. Offset(MyTotal,  0)  =  0 
MyTotal  =  MyTotal  +  1 

Loop  Until  MyFirstCell. Offset(MyTotal,  0).Address  =  MyLastCell. Offset(l,  0).Address 
Application. Calculation  =  xlAutomatic 

Dim  MyCount  As  Integer 
Dim  MyBestNode  As  Integer 
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Dim  MyBestValue  As  Double 


MyPrev.  Value  =  0 
MyNext.  Value  =  0 

Do 

MyBestNode  =  0 
MyBestValue  =  MyTarCell.  Value 

For  MyCount  =  1  To  MyTotal 

MyCurNode  =  MyCount 

If  ((MyFirstCell.Offset(MyCount  -  1,  0)  =  0)  And  (MyFirstCell.Offset(MyCount  -  1,-1) 
=  1))  Then 

MyFirstCell.Offset(MyCount  -  1,0)  =  1 
If  MyTarCell. Value  >  MyBestValue  Then 
MyBestNode  =  MyCount 
MyBestValue  =  MyTarCell.Value 
End  If 

MyFirstCell.Offset(MyCount  -  1,  0)  =  0 
End  If 

Next  MyCount 

MyFirstCell.Offset(MyBestNode  -  1,  0)  =  1 

MyPrev. Value  =  MyNext.Value 
MyNext.  Value  =  MyBestNode 

Loop  Until  MyBudget. Value  <  MyCost. Value 

MyFirstCell.Offset(MyBestNode  -  1,  0)  =  0 

End  Sub  'End  MyAveSolver 

VBA-MinCov  Algorithm 

Sub  MyMinSolver() 

Dim  MyCurNode  As  Range 

Set  MyCurNode  =  Worksheets("IO").Range("T31") 

Dim  MyFirstCell  As  Range 

Set  MyFirstCell  =  Worksheets("IO").Range("E2") 

Dim  MyLastCell  As  Range 

Set  MyLastCell  =  Worksheets("IO").Range("E674") 
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Dim  MyTarCell  As  Range 

Set  MyTarCell  =  Worksheets("IO").Range("T19") 


Dim  MyTarMinCell  As  Range 

Set  MyTarMinCell  =  Worksheets("IO").Range("T20") 

Dim  MyBudget  As  Range 

Set  MyBudget  =  Worksheets("IO").Range("N13") 

Dim  MyCost  As  Range 

Set  MyCost  =  Worksheets("IO").Range("N23") 

Dim  MyPrev  As  Range 

Set  MyPrev  =  Worksheets("IO").Range("T29") 

Dim  MyNext  As  Range 

Set  MyNext  =  Worksheets("IO").Range("T30") 

Dim  MyTotal  As  Integer 
MyTotal  =  0 

Application.Calculation  =  xlManual 
Do 

MyFirstCell.Offset(MyTotal,  0)  =  0 
MyTotal  =  MyTotal  +  1 

Loop  Until  MyFirstCell.Offset(MyTotal,  0). Address  =  MyLastCell.Offset(l,  0). Address 
Application.Calculation  =  xlAutomatic 

Dim  MyCount  As  Integer 
Dim  MyBestNode  As  Integer 
Dim  MyBestValue  As  Double 
Dim  MyBestMinNode  As  Integer 
Dim  MyBestMinValue  As  Double 


MyPrev.  Value  =  0 
MyNext.  Value  =  0 

Do 

MyBestNode  =  0 
MyBestMinNode  =  0 

MyBestValue  =  MyTarCell.Value 
MyBestMinValue  =  MyTarMinCell.Value 
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For  MyCount  =  1  To  MyTotal 
MyCurNode  =  MyCount 

If  ((MyFirstCell.Offset(MyCount  -  1,  0)  =  0)  And  (MyFirstCell.Offset(MyCount  -  1,-1) 
=  1))  Then 

MyFirstCell.Offset(MyCount  -  1,0)  =  1 

If  MyTarCell. Value  >  MyBestValue  Then 
MyBestNode  =  MyCount 
MyBestValue  =  MyTarCell.Value 
End  If 

If  MyTarMinCell. Value  >  MyBestMinValue  Then 
MyBestMinNode  =  MyCount 
MyBestMinValue  =  MyTarMinCell.  Value 
End  If 

MyFirstCell.Offset(MyCount  -  1,  0)  =  0 
End  If 

Next  MyCount 

If  MyBestMinValue  =  MyTarMinCell. Value  Then 
MyFirstCell.Offset(MyBestNode  -  1,  0)  =  1 
MyPrev. Value  =  MyNext.Value 
MyNext.Value  =  MyBestNode 
Else 

MyFirstCell.Offset(MyBestMinNode  -  1,  0)  =  1 
MyPrev. Value  =  MyNext.Value 
MyNext.Value  =  MyBestMinNode 
End  If 

Loop  Until  MyBudget. Value  <  MyCost. Value 

MyFirstCell.Offset(MyBestNode  -  1,  0)  =  0 
MyFirstCell.Offset(MyBestMinNode  -  1,  0)  =  0 

End  Sub  'End  MyMinSolver 
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Appendix  C 


The  solver  input  parameters  are  presented  below.  All  four  screen  which  allow  user  input 
or  selection  are  shown. 


Solver  Parameters  V7.1 


Set  Cell: 

Equal  To:  ♦  Max  Min 
By  Changing  Variable  Cells: 


Value  Of: 


$E$2:$E$674  <=  $D$2:$D$674 
$E$2:$E$674  =  binary 
$N$23  <=  $N$13 


$E$2:$E$674 

Subject  to  the  Constraints: 

Model 


Change 


Delete 


Solve 


Close 


Options 


|  OptQuest  Solver  - 

Add  Variables 


Reset  All 


Help 


_l 


Figure  25.  Solver  Objective,  Variable,  and  Constraint  Input  Screen 


Solver  Model 


Original  |  Transformed  j  Diagnosis  |  Options 


Unknown  Variables  Functions 

NonZeroes 

Close 

All  673  2 

Check  Model  | 

Smooth 

Check  For 

Gradients 

'  Structure 

C  Convexity 
(*  Automatic 

Quadratic 

Linear 

Bounds  673  Sparsity  % 

Integers  1 673  Total  Cells 

Solve  With  PSI  Interpreter  1  ♦  Excel  Interpreter 

Help 

Figure  26.  Solver  Interpreter  Selection  Screen 
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Figure  27.  General  Solver  Options  Selection  Screen 


OptQuest  Solver  Options 


General  Problem 


Max  Time: 


Iterations: 


1000000  seconds 


Load  Model... 


Save  Model., 


Precision  (Obj  Fun) :  j  0 , 00 1 
Precision  (Dec  Var):  [o.OOl 
Number  of  Solutions  to  Report:  [l 

F  Use  Same  Sequence  of  Random  Numbers  with  Seed  =  f 
[*  Check  for  Duplicated  Solutions 
I  Show  Iteration  Results 
[✓  Assume  Non-Negative 
F  Bypass  Solver  Reports 
r  Solve  Without  Integer  Constraints 

Stop  when  Objective  Hasn't  Improved  after  100  Iterations 


Figure  28.  OptQuest  Solver  Options  Input  Screen 
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